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Abstract 

^| In these lectures we discuss the Balitsky, Fadin, Kuraev, and Lipatov (BFKL) theory, 
which resums the leading logarithmic contributions to the radiative corrections to parton 
scattering in the high-energy limit, and we apply it to hadronic two-jet production at 
large rapidity intervals. 
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1 Introduction 



This is an expanded version of a few lectures I gave at the Universita' di Salerno and 
Torino, Italy, on the perturbative QCD pomeron and on jet physics at large rapidity 
intervals. The outline of the lectures is as follows: in the Introduction the parton- 
model and the factorization pictures, propedeutic to any calculation of strong-interaction 
processes, are sketched; in sect. 2 two-jet production at hadron colliders at the leading 
order in a s is discussed, by examining the parton kinematics and dynamics in the exact 
configurations and in the large-rapidity limit; in sect. 3 higher-order corrections to two- 
parton scattering are considered in the limit of a strong rapidity ordering of the final-state 
partons. In resumming the leading logarithmic contributions to the radiative corrections, 
the BFKL equation describing the gluon-ladder evolution in transverse momentum is 
introduced; in sect. 4 the BFKL formalism is applied to the description of inclusive 
two-jet production at large rapidity intervals. 

1.1 Jet Production and Factorization 

The importance of jets in hadronic collisions has been realized since the conception of the 
parton model, where jets arise from the scattering between the constituent partons of the 
colliding hadrons. On the experimental side, hadronic jets were then actively searched 
at the CERN ISR collider, but their existence was demonstrated only at the CERN SPS 
collider |IJ, by applying a cut in the total scalar transverse energy of the event. This 
strongly reduced the soft-hadron background of low transverse-energy partons due to 
the "underlying event", i.e. to the soft scattering between the spectator partons, and 
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left mainly high transverse-energy jets [] 

In the parton model the scattering process is factorized into two regions: i) a short- 
distance region, characterized by a hard scale Q, of the order of the jet transverse energy 
~ 10 GeV and thus typical distances of about 10~ 2 fm, which describes the primary scat- 
tering between the partons; ii) a long-distance region, characterized by a hadronization 
scale ~ Aqcd and distances of about 1 fm, which describes how the scattering partons 
split from the parent hadrons and how the scattered partons hadronize. 

To see intuitively how this factorization comes about in the parton-model picture, let 
us imagine to sit in the center-of-mass frame of the hadron scattering. Then we observe 
two kinematic effects: i) each hadron looks like a pancake, Lorentz contracted in the 
direction of the collision. Accordingly, the time it takes to parton a(b), within hadron 
A(B), to go through hadron B(A) will be Lorentz contracted; ii) the parton interactions 
within each hadron are time dilated. Therefore, the time it takes to a scattering parton to 
go through the other hadron is much less than the interaction time between two partons 
in the same hadron. That is, parton a (b) sees hadron B(A) as a frozen distribution of 
partons. Besides, the momentum transfer Q being large, parton a(b) probes on hadron 
B(A) a region of transverse size ~ within which, unless the hadron is densely 

packed, it will meet only parton b(a). 

Then the production process can be thought of as the convolution of two contributions 

V = Yjf dx A J dx B fa/A{xA^ 2 )fb/B{xB^ 2 )6 a b{(A/xAi(B/xB,as{^ 2 ),Q 2 / ' ^) , (1) 
ab ^ A ^ B 

where a and b are summed over the parton species and flavors; (a, (b and a are 
a An introduction to the experimental discovery and to the theory of hadronic jets is given in ref.g]. 
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parametrized by the kinematic variables of the produced hadrons or jets; a is the prob- 
ability for the hard scattering to happen, which can be computed through perturbative 
quantum chromodynamics (QCD) calculations; f a /A{ x A, H 2 ) is the probability that par- 
ton a(b) is found in hadron A(B) carrying a fraction xa(b) of the parent hadron mo- 
mentum. fa/A^XA, H 2 ) is a universal quantity, i.e. it does not depend on the particular 
hard scattering a considered, it can not be computed in perturbation theory and must 
be derived from the experimental data. 

Eq.(0) expresses the factorization in the parton model, which is valid in the region 
where x is fixed and Q 2 — > oo. For finite, and large, values of Q, the production cross 
section is expanded in powers of l/Q; then eq. ([!]) may be shown to be valid for the leading 
term of the expansion ||, called the leading twist from operator product expansion 
terminology. The remainder of the expansion, which does not usually factorize, is called 
the higher twist. 

Then we proceed to calculate the partonic cross section a, as an expansion in the 
strong coupling constant a s 



a ((/x,a s ,Q 2 /iJ 2 ) = a 



..h 



71=1 



(2) 



where the power h and the coefficients depend on the production process. In doing 
the calculation, we set the parton masses and transverse momenta to zero, since their 
contribution to eq.(0), of 0(Aq C d/Q), is in the higher twist. Calculating the coefficients 
with i > 1, we encounter both ultraviolet and infrared divergences. The ultraviolet 
divergences are due to the virtual radiative corrections, they are an artifact of the pertur- 
bative expansion and are subtracted away by using a renormalization prescription. This 



introduces in the calculation the dependence on a renormalization scale /i. Left over are 
then the infrared divergences, divided into soft divergences, which appear when a parton 
momentum vanishes, and collinear divergences, due to the collinear production of mass- 
less partons in a vertex. After including all the virtual and real radiative corrections to a 
given coefficient the soft divergences cancel out. The left-over collinear divergences 
are universal, i.e. are not specific to the coefficient o-w we obtained them from. They 
are the outcome of the collinear emissions from parton a(b) on its evolution toward the 
hard scattering. This evolution is space-like, hence during the evolution the absolute 
value of the parton virtuality grows. In particular the collinear emissions describe the 
distribution of parton c(d) of momentum fraction Ca(B) and virtuality /i 2 , within parton 
a(b) of momentum fraction xmb) an d zero mass and transverse momentum. Since the 
collinear parton evolution is universal we factorize it into the parton density f a /A{xA, /-f 2 )- 

While both the parton density and the partonic cross section depend on the renormal- 
ization/f actor ization scaleQ /i, their convolution (|]) does not, since the physical process 
does not depend on the value of the virtuality fi 2 we stop the collinear parton evolution 

at , 

= < 3 > 

Replacing into eq.(|3]) the factorization formula (|J), we obtain renormalization group 
equations for the coefficients and for the parton densities ||. The ones for the 
parton densities 



''Even though they are often chosen to be the same, the renormalization and factorization scales are 
in principle unrelated (see sect. 14.3 of ref.fl). 



are called Altarelli-Parisi, or DGLAP, evolution equations ||, and describe the distribu- 
tion of parton c within parton a, as pictured in the paragraph above[]. The functions P ac 
are known as Altarelli-Parisi splitting functions and may be computed in perturbative 
QCD as an expansion in a s 

^ C = Eg) p^\0- (5) 

Eq.(f|) resums the collinear logarithms in the evolution at a given accuracy, determined 
by the splitting functions (|5|). Namely, the leading order (LO), or one-loop, splitting 
functions || yield the (all-order) leading logarithmic evolution in eq.([|), the next-to- 
leading order (NLO), or two-loop, splitting functions [0] yield the next-to- leading loga- 
rithmic evolution in eq.(|4]), and so on. Accordingly, the running of a s in eq.([|) must be 
determined to the corresponding loop-accuracy. 

Eq.(|D states the indipendence of the physical cross section on the factorization scale 
/1 2 , and it holds exactly if we know the full expansions @) and (||). However, in a fixed- 
order expansion of the factorization formula (|l|), say at eq.(|3|) holds only up to 
corrections of 0(a™ +1 ). Hence the more terms we know in the expansion of the partonic 
cross sections (0) and of the splitting functions the less our evaluation of eq.([I|) 
depends on the unphysical scale fi 2 . 

In inclusive jet production the partonic cross sections are known at best at NLO 
in a s . Namely, at present 4-parton || and 5-parton |J NLO matrix elements are avail- 
able. By combining the 4-parton NLO matrix elements [BJ with the 5-parton LO matrix 



elements ||10|| , NLO one-jet [Ly] and two-jet [ 12|— [ 14| inclusive distributions have been 



c An introduction to eq. (Q) and its relation to the renormalization group may be got from the horse's 
mouth |q|. 



computed. Besides reducing the dependence on the factorization scale /i 2 , they allow 
one an analysis of the dependence of jet production on the jet-cone size and of the par- 
ton distribution within the cone Jll| . They appear to be in good agreement with the data 
on one-jet |T3] and two-jet [16|]-|[l!J inclusive distributions from the Fermilab Tevatron 



Collider. 

1.2 Large Rapidities and the BFKL Evolution 

In section |1.1| we have outlined how the collinear factorization applies to an inclusive 
production process and we have quoted the main results of the standard hadronic jet 
analysis. We have assumed that in the production process there is only one hard scale Q, 
the momentum transfer, of the order of the jet transverse energy; and that the hadron 
center-of-mass energy yfs is of the same order as Q. However, at the Fermilab Tevatron 
Collider and at future hadron colliders the semihard region of the kinematic phase space, 
where \fs ^> Q, is accessible. If \fl = \/xa%bs is the parton center-of-mass energy, then 

The logarithms, ln(l/x), appear in the evolution of the parton densities; the logarithm, 
ln(s/Q 2 ), parametrizes the hard scattering a. In the semihard region the left-hand side 
of eq.(||) is large. Then in the production process we may have large logarithms, ln(l/x), 
or large logarithms, ln(s/Q 2 ), related as we will see to the occurrence of large rapidity 
intervals between the produced jets. Thus, large non-collinear logarithms may have to 
be resummed either in the splitting functions (|J) or in partonic cross section (|J). A 
unified treatment of factorization which takes into account the general case where both 
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these logarithms are large, and at the same time includes the collinear factorization ([!]), 



is given in ref. |20 



In these lectures, we deal with jet production in the semihard region, at large momen- 
tum fractions x's of the incoming partons In this case there are no large logarithms, 
ln(l/x), in the splitting functions ([5]). Thus the parton densities evolve according to the 
usual Altarelli-Parisi evolution and the collinear factorization (|I]) is suitable to describe 
the production process. However large logarithms, \n(s/Q 2 ), may appear in the partonic 
cross section (|2|), and we consider the problem of resumming them. 

Resummation techniques of large logarithms, ln(s/Q 2 ), date back to the studies of the 
high-energy limit of quantum electrodynamics (QED) ||22|| . It is easy to see using power- 
counting arguments that at high center-of-mass energy yfs and fixed momentum transfer 
Q, the leading contribution to a given scattering amplitude comes from photon exchange 
in the crossed channel t. This contribution is responsible for letting the total cross section 
approach constant values at very high energies, in accordance with Pomeranchuk theorem 
|[23|| . Let us suppose then that along the photon exchanged in the t channel a fermion 
pair is emitted. The ensuing 0(a 2 ) radiative corrections to the scattering amplitude 
may contain large logarithms, ln(s/Q 2 ) P2J. The leading powers of [a 2 \n(s/Q 2 )} may 
be resummed in the perturbative expansion, and the resummed series yields the total 
cross section for the amplitude with exchange of one photon in the t channel. 

The same techniques then have been applied by Lipatov and collaborators to non- 
abelian gauge theories [pl}| -[P6|, and in particular to perturbative QCD p7| . Like in 
QED, in the limit of high center-of-mass energy \fl and fixed momentum transfer Q 



7 



the leading contribution to a given scattering amplitude in a physical gauge comes from 
gluon exchange in the t channel. However, due to the gluon self-coupling interaction, 
the radiative corrections which contain large logarithms, \n(s/Q 2 ), appear at 0(a s ). 
Therefore, they are roughly a factor a s /a 2 stronger than the corresponding corrections 
in QED, and may be of practical importance in computations of high-energy scattering. 

2 Two-jet production at leading order 

In this section two-jet production at hadron colliders at the leading order in a s is dis- 
cussed. We begin by defining the jet kinematic variables by which the high-energy 
process will be described; the two-parton kinematics as a function of the jet variables 
and the factorization formula that connects the parton subprocess to jet production are 
introduced; the parton dynamics is discussed and the relevant subprocess in the large- 
rapidity limit is computed; two-jet production at LO at the Tevatron and the LHC 
colliders is calculated using the kinematics and dynamics in the exact form and in the 
large-?/ approximation; finally, it is sketched how two-jet production can be used to gain 
knowledge about the parton densities. 

2.1 Rapidity 

Let us consider two inertial frames with the longitudinal axes oriented along the beam 
direction of the colliding particles, and in relative motion in the beam direction. Then 
the motion of a particle, at rest in a frame, is characterized in the other frame by the 



8 



relation 



tanhy = ^, (7) 



where E, p and y are respectively the particle energy, longitudinal momentum and 
rapidity Inverting eq. ([?]), we obtain the definition of rapidity as 

1 E + Pu 

y=t\ n ^X, ( 8 ) 

Eq.(^) states that energy and longitudinal momentum scale like coshy and sinhy, respec- 
tively. The constant of proportionality is given by the mass-shell condition. Introducing 



the transverse mass, m± = yjm 2 + p\, where m is the particle mass and p± the absolute 
value of the momentum in the plane transverse to the beam direction, the mass-shell 
condition reads E 2 = p 2 + m\ . Using it in eq. (^) , the particle energy and longitudinal 
momentum are expressed by 

Vu — rn±smh.y] (9) 
E = m±coshy. 

It is often convenient to introduce light-cone coordinates, by combining energy and 
longitudinal momentum 

P ± = E±p l] , (10) 

then the scalar product of two vectors is p ■ q = (p + q~ + p~q + )/2 — p± ■ q± and the 
non- vanishing components of the metric tensor are 

2#+- = 2#_+ = -g xx = -g yy = 1. (11) 

Using eq . (|To|) , the particle 4- momentum is parametrized by mass, transverse momentum 
and rapidity as 

p = (m ± e y , m±e~ y ; p ± ) , (12) 



where = (p± cos<p,p± sin0) and <fi is the azimuthal angle between the vector p± and 
an arbitrary vector in the transverse plane. For massless particles we have from eq.(0) 
that tanhy = cos 6, with 9 the angle between the directions of the scattered particle and 
the beam. Then eq.([8|) transforms to 

i , i + cos e , e , , 

y(m = 0) = 2 ln T^o^ = ~ lntan 2- (13) 

Eq.fll3l) defines the pseudo-rapidity, r\ = y(m = 0). This, and not the rapidity y, is the 
variable used in high-energy experiments, since its determination requires just tracking 
the particle direction. For massless particles, or in the case where p± ^> m, the two 
definitions coincide. However, it is worth reminding the difference, since it is the rapidity 
y to define the correct measure of the particle phase space 

d 3 p dy d 2 p ± 
2£(2tt) 3 47r(27r) 2 ' 1 ' 

and to transform additively under boosts in the beam direction. Thus, the shape of a 
particle multiplicity distribution dN/dy is boost invariant. To see the consequences of 
considering, on the contrary, multiplicity spectra in pseudo-rapidity dN/drj, let us take 
a typical hadron-hadron collision event, where the low transverse-momentum particles 
from the underlying event fill up, rather uniformly, the phase space in y. Now, consider 
eq.(||), and the analogous relation for pseudo-rapidity sinh^ = p ./p±. Massive particles 
with p± m will have \r]\ > \y\, hence they are pushed away from the pseudo-rapidity 
central region t] ~ 0. Then a fake excess of massless particles (photons) at rj ~ is found 
in the pseudo-rapidity distribution dN/ dr] ^8] . 
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2.2 Elastic Scattering 



Let us consider the scattering of parton a in hadron A with parton b in hadron B, in 
the hadron center-of-mass frame, which we choose to be the lab frame. In light-cone 
coordinates, the 4-momenta of the incoming partons are 

Pa = (Vsza,0;O), (15) 
Pb = (o,y/sx B ]0). 



The momenta of the outgoing partons, which at LO we identify with the observed jets, 
are given according to the parametrization ( |T2D by 

k=(k i± ev%k i± e- y <;k i± ), (16) 

with i = a,b. We introduce then the rapidity boost y = (y a +yb)/2, i.e. the rapidity of the 
parton center-of-mass frame with respect to the lab frame, and the rapidity difference 
y = 2y* = y a — yt,, which is twice the rapidity y* of a parton in the center-of-mass 
frame. Since in jet production the parton final states are indistinguishable, without 
losing generality we assume that y > 0. 

Transverse momentum conservation requires that k a ± = kb± = k±, while light-cone 
momentum conservation determines the momentum fractions x's of the incoming partons 
as a function of the jets kinematic variables 

x A = %e Va + %e y » = 2^eV C oshy*, 
V s V s V s 

(17) 

x B = %e~ Va + -4=e- m = 2^=e-y C oshy*. 
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Taking the ratio, we obtain 2y = \b.(xa/xb), and in agreement with the definition of 
rapidity boost, y = when the center-of-mass frame is at rest in the lab frame. The 
Mandelstam invariants, which characterize the parton scattering, are given in terms of 
the jet variables by 

s = Ak\_ cosh 2 y*, 

t = -2k 2 ± coshy*e- y *, (18) 
u = —2k 2 L coshy*e y . 

Note that the Mandelstam invariants do not depend on y, since y does not belong to the 
center-of-mass frame kinematics. Using in eq. fll8|) the relation between the Mandelstam 
invariants and the scattering angle 9 in the center-of-mass frame, 




u = --(1 + COS0), 

we find y* = — lntan(#/2), in agreement with the definition ([13|) for massless partons. 
2.3 Two-jet production 

Having described the parton kinematics, we have to establish how the parton scattering is 
related to jet production at the hadron level. To do that, we need a factorization formula, 
as outlined in the introduction. The factorization formula for two-jet production cross 
section in terms of the jet rapidities and transverse momenta is in general given by 

d 2 k a ±d 2 k b _Ldy a dy b ? / dxAXB f*/ A ( XA > ^ )Ij/b{^b^ ) ^^^2^^^' ( 20 ) 
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where i and j are summed over the parton species and flavors. The parton cross section 
for elastic scattering is 

(2n) 4 5 4 (p a +p b -k a - k h ) dy a d 2 k a± dy b d 2 k b± 2 
^ = 2s M^F ' ijl ' ( } 

where we have used the parametrization ([[4]) for the phase space of the outgoing partons, 
and where |Mjj| 2 are the squared matrix elements, summed (averaged) over final (initial) 
parton helicities and colors^. 

Using eq.(|T7p in the momentum-conserving 5-function to fix the parton momentum 
fractions, and replacing eq.([H]) into eq . (|20D , we obtain the factorization formula for 



two-jet production at LO 



J2 x Afi/A(x A , ii 2 ) x B f j/B (x B ^ 2 )^-, (22) 



dk]_dy a dy b ij dt 
where the cross section for elastic scattering is 



do a _ \Mij 



|2 



(23) 



dt 16ns 2 

As noted after eq. (|l8D , the Mandelstam invariants, and so |M i:) '| 2 , depend on the trans- 
verse momentum kj_ and the rapidity difference y, but not on the boost y\ while the 
parton momentum fractions x's depend also on y (cf. eq. (pp). This affects the choice 
of variables which characterize two-jet production. To help visualize this, let us consider 
in Fig. m the plane of the jet rapidities. Then if the object of our study is the parton 

scattering, it is convenient to either fix y (the solid lines of Fig. [I]) or integrate it out, 

d In eq.([2l|) we have not included the symmetry factor 1/2 for identical final-state partons. We will 
keep that into account in the phase space with y > by summing only over the momentum configurations 
of different final-state partons. This introduces no doublecounting, since everything is symmetric under 
the exchange of the two jets. 
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Figure 1: Plane of the jet rapidities. The solid curves are lines of constant y, the dashed 
curves are lines of constant y. 

since its variation induces just a varying contribution from the parton densities. If in- 
stead we are interested in the parton densities, then it is convenient to fix k± and y (the 
dashed lines of Fig. |l|), thereby fixing the contribution of the parton scattering to two-jet 
production. 



2.4 Parton Dynamics at Large Rapidities 

In order to use the formula for two-jet production at LO (f22|), we need to compute the 
squared matrix elements |Mjj| 2 . They have been evaluated in ref. [29], and are 6 



2 4 s 2 + u 2 



= \M q ?-. q ?\ 2 = g— , (24) 

±(s 2 + u 2 s 2 + t 2 \ 8 s 2 
l^—l 2 - g — + "27 S . (25) 



e The explicit calculation for the subprocess qg — > qg can be found in sect. 7.2 of ref. [p0[ 
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\-M.qg-*g' q'\ 
\M- q q-+gg\ 

\M-gg^qq\ 
\Mqg^qg\ 



\M 



9 9^9 9 \ 



A ( 5 2 _i_ f. 


^ f 2 + ? "} 2 \ 


Q \ + 2 


S / 


At -\- u 




9 s 2 ' 




32t 2 + « 2 


8 t 2 + M 2 


27 


3 s 2 ' 


1 f 2 + u 2 


3t 2 + u 2 


6 £m 


8 s 2 ' 


s 2 + u 2 


4s 2 + « 2 


i 2 


9 su ' 


9 / £u 




2 3 "F 





u 



27 ; 



(26) 

(27) 
(28) 
(29) 
(30) 
(31) 



where q and q' are quarks of different flavors, and |Mjj| 2 = \M.ij\ 2 , with g 2 = 4iTa s . 



Now, we examine the dependence of the parton subprocesses ( P4] - |3T|) on the rapidity 
interval y = 2y*, and as a sample we consider the subprocess g g — > g g ([H]). Replacing 
in it the Mandelstam invariants (IH), we obtain 



\M 



9 (4cosh z y*-l 
"2 4 cosh V 
then the corresponding parton cross section (|23"D is 



(32) 



do", 



99^99 



9 2 1 (4 cosh 2 - 1 



- 7r a 



dt 2 b k ± ^4 C osh 2 ?/*) 

At small rapidities, eq. (|33|) increases with y, 

at 3 



(33) 



(34) 



Then we consider the large-?/ limit. For jy ^> 1, we find from eq.flHf) that 



(35) 
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Combining eq.(55|), the rapidity difference y at large rapidities becomes 

2/ = ln(-0, (36) 

i.e. the large-y limit is equivalent to the high-energy limit at fixed momentum transfer. 
Four Feynman diagrams contribute to the subprocess g g — > g g: the exchange of a gluon 
in the t channel (FigfJ); the exchange of a gluon in the s(u) channel, obtained from Fig|| 
by crossing leg k a with leg Pb(kt); the four-gluon coupling. Since s ~ |w| ^> \t\, only the 
diagram with exchange of a gluon in the t channel contributes in the limit y ^> 1, as 
long as we consider only physical polarizations for the external gluons. 



Pa ^a a 1 


10 ka -^a' a ' 




*° q, c 


p b A b b | 





Figure 2: g g — > g g scattering, with exchange of a gluon in the i channel. We label each 
external line with momentum, helicity and color. 

Inverting the metric tensor ( |TTD in order to let it act on the helicity labels of a 
complete basis of (polarization-like) unit vectors 

<T = v»g xx 'v» y = 2«t£ + - < • vl, (37) 

and using eq. (|T5|) in the light-cone vectors 

v + = (1,0,0) = — ^p a , (38) 
x A \/ s 

= (0,1,0) = —^-=p b , 

x B V s 
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we can rewrite the metric tensor as 

gfl u = 2 PtPb+Patf _ ^ ( 39) 

where 8^ is a Kronecker delta over the transverse components. Next, we compute the 
amplitude for the exchange of a gluon in the t channel, to leading 0(s/\t\), and use 
eq . (|39D in the gluon propagator. Then the amplitude of Fig.[| becomes 



^Cv'vv - 9,f°°" [w,(f. + t)» + V..(-». + ?k-9,...(4 + P.)" 
v ' s t 

■ 9s f Wc [g^ b , {p b + h) P - sw {h + g) w + g„ bP (q - p b Y b '] (40) 
• ^iPa)e^{p b )e^{K)eZ{hl 

~ -2iglr' c g» a »/ 1 f Wc g w e^{ Pa )e^ 

where f abc are the SU(N C ) structure constants, with N c = 3 the number of colors, and 
are the gluon polarization vectors of helicity A, and q 2 = t. Since in eq.([K]) we are 
considering physical polarizations, e\(p) ■ p = 0, the leading contribution in s/\t\ comes 
from combining the helicity-conserving terms in the 3-gluon vertices with one of the two 
light-cone polarization modes in the gluon propagator. Then, we square the amplitude 
and sum over helicities and colors. We replace the sum over the gluon helicity states by 

^ \ n ■ p [n ■ p) 1 J 

where n is an arbitrary 4-vector. This is equivalent to use an axial gauge. For example, 
a convenient choice for the sum over the helicity states of gluon p a is to take n = p^. 
Then eq.(^T|) becomes 

E CM) tZiPa) = - (g" aUa - p " a pla + p " rt a \ = sr a , (42) 

X a \ Pa'Pb J 
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where the use of eq.(^) shows that only physical polarization states are summed. Using 
eq. fl42|) in the amplitude (fK]) squared, we obtain for the sum over the helicity states of 
gluons p a and k a , 



Kt 



2 l + O I- , (43) 



which shows that helicity is conserved at the jet-production vertices, up to subleading 
terms. 



The procedure we outlined in eq. (|40| - f43|) is not gauge invariant, since the terms we 



neglected in eq. (|40D are subleading only if we consider physical polarizations, that is if 
we use a physical gauge. In order for gauge invariance to hold, Ward identities must 
be fulfilled, i.e. by replacing in the scattering amplitude one or more of the physical 
polarizations with the longitudinal one, e^(p) — > p^, the amplitude must vanish^]. It 
is easy to see that this is not the case in eq.(^). To obtain then a gauge- invariant 
amplitude, one must consider the full expression of the amplitude for the exchange of 
a gluon in the t channel, given in eq. (^0|) , and include also the Feynman diagrams for 
gluon exchange in the s or u channel and for four-gluon coupling. Thus, amplitude ( |40|) 



assumes the form p5 



-2< 9t f aa c r Ma , al j r c r w w e£ ( Pa ) ( Pb ) e£ (k a ) (k b ) , (44) 



with 



Pa-Pfe 2(p a -p b ) 21 



(45) 



■'For a detailed discussion of gauge invariance and Ward identities in non-abelian gauge theories, see 
sect. 8.5 ofref.il, or sect. 7.2 of ref.EJ. 
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and the analogous expression for r M& Mi/ , by exchanging the labels a and b. It is straight- 
forward to check that amplitude ([14]) is gauge invariant, up to subleading terms. Besides, 
we could replace the sum over the gluon helicity states by 



-<r, 



(46) 



and obtain the right helicity counting, 



Aa 



(47) 



However, considering only diagrams with gluon exchange in the t channel, as we did in 
eq. (f40D , is a correct procedure as long as we work consistently in a physical gauge. 



The sum over colors then yields 



jaa'c ja'ac' jbb'c jb'bc' ^^y2 



(48) 



where Ca = N c is the Casimir factor of the adjoint representation of SU (3), under which 
the gluons transform. Replacing eq.([|3|) and ([i"8"D in the amplitude (|40|) squared, and 
averaging over the initial helicity and color states, we obtain 



IM, 



2 4Cj 4 s 2 9 4 s 2 



(49) 



which is in agreement with eq . (|3l"l) , if there we take the Mandelstam invariants in the 
large-rapidity limit (|35|). 



Then we examine the other subprocesses that contribute to two-jet production. Re- 
placing eq.(p5[) into eq. ( p^ - pO| ) , we realize that subprocesses ([27]), ( |28| ) and (^) do not 
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give a leading contribution in (s/|t|), since they have no diagrams with gluon exchange 
in the t channel. Subprocesses fl24]) , (j^), (|26| ) and (|30D yield, in the large-rapidity limit, 



s 2 

,4 6 



|M„,^„^ = |M g( ™| 2 = \M qq ^ qq \ 2 = -g^, (50) 



s 2 



|M gs _> g ^ = 2g\ j2 . (51) 

Eq. (^Dj) and (^T|) may be obtained by taking in each subprocess the leading contribution 
to the diagram with gluon exchange in the t channel. Taking the ratios of eq . (|49|) , (]50l). 
and (|51|), we obtain 

\M gg ^ 9g \ 2 = - A \M qg ^ qg \ 2 = {-\ 2 \M qq ^ qq \\ (52) 



A 

where 9/4 = Ca/Cf is the relative color strength in the jet-production vertices (CV = 
(N 2 — 1)/2N C = 4/3 is the Casimir factor of the fundamental representation of SU(3), 
under which the quarks transform). Thus it is enough to compute one of the three sub- 
processes in eq . (|52|) to describe the parton dynamics in the large-rapidity limit. Having 
done it for the subprocess gg — > gg, we can include the contribution of the others in 
two-jet production at LO ([22"D by using the effective parton density [[31], 

f eff (x, /i 2 ) = G(x, /i 2 ) + | E [<?/(*. M 2 ) + Qf(x, fi 2 )} , (53) 

where the sum is over the quark flavors. We dropped the hadron label in eq. (|53|) , since 
the effective parton density is charge-conjugation invariant. To complete the discussion 
of the parton dynamics at large rapidities, we note by using eq.(|23|) and (|4"9"D, or taking 
the large- y limit of eq. (|3"3"D , that the parton cross section does not fall off as the parton 
center-of-mass energy rises, 

da 9 2 1 . . 

— = -jca 2 — . (54) 

dt 2 s t 2 K J 
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This is in agreement with expectations of a constancy of the total cross section at high 
energies 
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Also the kinematics simplify in the large-?/ limit: the parton momentum fractions 
([!]]) become 

x° A = ^e* (55) 

r ° - k± p-Vb 
x B — _e , 

i.e. each parton momentum fraction is determined by the kinematic variables of one jet 
only. Finally, we examine how the factorization formula (|22|) tranforms in the large-?/ 
limit. The phase space for the production of two partons in eq.(^l|) is 

f dy a d 2 k a± dy b d 2 k b _L 4 4 
^2 = / , /„ N 2 , /„ N 2 (2vr) 5 (p a +p b -k a - k b ) . (56) 



47r(27r) 2 47r(27r) 2 

We can fix the rapidities using the light-cone momentum conservation ([T7|) and the 
parametrizations (|T5| ) and (|16"D, 



2 y dy a dy b 5 (<fsx A - fc^e"" - fc^e^) 5 (v^s - &j_e Va - k ± e~ 
2 cosh (y a - y b ) + 1 



-;</(, 



(57) 



s sinh (y a - j/ 6 ) 

where the overall factor 2 on the left-hand side comes from the Jacobian in the change 
of variables (fLOf). Taking the large-;?/ limit in the right-hand side or using the light-cone 
momentum conservation at large y (|55|) in the left-hand side, eq.flSTD reduces to 2/s. So 
the phase space for the production of two partons (|56"D becomes, in the large-y limit, 

p -/5WM w!i,fo+w ' (58) 
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and the cross section for elastic scattering ( pT|) reduces to 



daij _ \Mij 



|2 



(59) 



which coincides with eq . (^3|) . Using then eq.(^) the factorization formula (|22|) becomes 



5^ = 4/^)4/^") f«. (») 

Note that we could have obtained the factorization formula ( |60|) from eq.(22) by simply 



recalling (cf. eq . (|35|) ) that in the large-?/ limit t — > —k\. However, the procedure we 
have followed above is suitable for generalizing eq.(|60"D to higher orders in a s . 



2.5 Phenomenology of two-jet production at leading order 

We are now ready to perform a quantitative study of two-jet production, as a function 
of the rapidity interval y between the jets. First, we mention that the dominance of t- 



channel gluon exchange as y increases, as described in sect. 2.4, has been verified through 



a calculation of two-jet production as a function of y at LO and NLO in a s [12], which has 



been found in good agreement with the corresponding data from the Tevatron collider 
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In this section we outline how to perform a LO calculation of two-jet production, 
collecting the information we have acquired in the previous sections, and we examine the 
accuracy of the large-rapidity approximation at LO, at the Tevatron proton-antiproton 
(y 7 ? = 1.8 TeV) and at the LHC proton-proton (^/i = 14 TeV) colliders. For an exact 
calculation we need the factorization formula (0), with the parton cross section given by 
eq.(|23D and ( p4"| - |31| ). The Mandelstam invariants and the parton momentum fractions are 
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given in eq.flTTD and (JTM) as functions of the jet kinematic variables. For the calculation 
in the large-?/ approximation, the factorization formula (|60"D is to be used, with the 
effective parton density (|53|), and the parton cross section given by eq. (|23|) and (|49|). 



The corresponding Mandelstam invariants and parton momentum fractions are given in 
eq.(H) and (H). 

Since we are mainly interested in understanding how the large-?/ approximation works 



on the parton dynamics, we fix the rapidity boost y = (cf. sect. |2.2| and [2.3|) . Then 
we take p± > 20 GeV, so that the parton densities are evaluated at x <; 2 ■ 10 -2 at the 
Tevatron (3 • 10~ 3 at the LHC). The choice of the factorization scale /i 2 is arbitrary (cf. 
sect. [TTTD ; we fix its value at the jet transverse energy, /i 2 = p 2 ± . For the LO evolution 
of the parton densities with x and /i 2 , we take the CTEQ parton densities f3"2fl . Finally, 
we evaluate a s (/i 2 ), scaled from a s (m^) = 0.12, using the one-loop evolution with five 
flavors. 

In Fig. |] we plot the two-jet production cross section da/dydy as a function of y, 
at the energies of the Tevatron and LHC colliders. The dashed (dot-dashed) curves are 
obtained from the exact (large-zy) evaluation of the cross section. As expected the large- y 
approximation improves as the size of the rapidity interval grows; at y ^ 3 the difference 
between the curves is about 20%, thus, it is not bigger than the theoretical uncertainty 
due to the choice of factorization scale. The dot-dashed curves of Fig. |] (the dashed ones, 
for large y) fall off with y. This is due to the parton densities. Indeed at fixed y the x's 
grow linearly with the jet transverse momentum (cf. eq. (|17|) and (|55|) ), thus at larger 
values of y the integration over p± starts at larger values of the x's, where the parton 
densities give a lesser contribution |33| , |34j] , [^5| . The behavior of the curves at small y 
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is different, though. From eq.flSl) we know that the parton cross section in the large-?/ 



approximation is constant as y grows. This, together with the reduction in phase space 
from the kinematics, yields the monotonic fall-off of the dot-dashed curves. The exact 
parton cross section fl33|), instead, grows with y at small rapidities (|34l). This yields the 
initial rise of the dashed curves, followed then by the parton-density suppression. 

100 



10 



.01 



p ± > 20 GeV 



s = 14 XeV-^ 



Vs = 1.8 TcV 



.001 

2 4 6 8 10 12 

J 

Figure 3: Two-jet production at the Tevatron and LHC colliders, as a function of the 
rapidity interval y, at y — 0. The dashed and dot-dashed lines are respectively the exact 
and large-?/ LO cross sections. 

Next, we sketch how two-jet production may be used to obtain information on the 
parton densities. For x ^ 10~ 2 the quark densities are well known from precision mea- 
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surements of the structure function F2 in deeply inelastic scattering (DIS) |56]| ; from 
the target isospin dependence in Drell-Yan scattering ||37|| ; from the W ± asymmetry in 



hadron collisions p8[ . The gluon density, though, is known only at x ~ 0.4 ]39]]. For 
x < 1CT 2 the gluon density is presently extracted from the quark sea in measurements of 
F2 in DIS [ 40| . Also two-jet production may be used to determine the gluon density [14}|, 



|19| . This can be achieved by considering two-jet production with the jets on opposite 
sides (OS) on a plane in rapidity and azimuthal angle, the lego plot, or with the jets on 
the same side (SS) of a lego plot. 

In OS two-jet production we fix y = and let y = 2y* grow; the center-of-mass 
frame is at rest in the lab frame, and the incoming partons carry the same fraction of 
the parent hadron energy, xa = xb- From eq.(^), the cross section behaves like 

<?OS ~ J2 fi/A(%A, /i 2 ) fj/B^A, ^) &OS- (61) 

At the Tevatron collider xa is typically ^ 10~ 2 , and its lower bound grows with y, as 
described in the paragraph above. Since at these values of xa the quark density is well 
known [36|-|j38|, eq.(|6lj) may then be used to extract information on the gluon density 



at the same values of xa- 

In SS two-jet production we fix y = and let y = y* grow; then one x, say xa, 
grows and the other, xb, decreases, i.e. one of the incoming partons is rather energetic 
and the other is wee, and the center-of-mass frame is boosted by y in the direction of 
the energetic parton. The wee parton is predominantly a gluon, since the gluon density 
dominates at small x. Then the cross section behaves like 

oss ~ U/A&A, ^ 2 ) g(xB, /U 2 ) a ss , (62) 
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with as small as 2 • at the Tevatron collider. Since xa is now very large, the 
parton density fi/A(xA, H 2 ) is dominated by the valence quarks, and its value is well 
known. So eq.(|62|) may be used to gain information on the gluon density at small values 
of x. 



3 The perturbative QCD pomeron 

In this section higher-order corrections to two-parton production are considered, in the 
limit of a strong rapidity ordering of the final-state partons. This kinematic regime yields 
the leading logarithmic contributions, of the type ln(s/|t|), to the radiative corrections. 
Tree-level three-gluon amplitudes, and the related production cross section, and one-loop 
elastic amplitudes are computed in this limit. Multigluon amplitudes are then introduced 
and used to evaluate the elastic scattering amplitude with exchange of a two-gluon ladder 
in the t channel. In doing this, the BFKL equation which describes the gluon-ladder 
evolution in transverse momentum is derived^. The solution is found for color-octet |25[ . 



and for color-singlet exchange |26|. The singlet solution is related through unitarity to 



the total parton cross section, with exchange of a one-gluon ladder in the t channel [p?J . 



3.1 The li-parton kinematics 

Let us assume that in the parton scattering n + 2 partons are produced (Fig. ||). The 

4-momenta of the incoming and outgoing partons are parametrized like in eq . fll5|) and 

s The BFKL equation also describes the small- a; evolution of the gluon density, and in that context 
it has been derived in ref . Q , Q . 
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(|16|), with % = 0, n + 1. Momentum conservation requires that 



n+l 



o = £k i± , 



i=0 

n+l u 

ft?. 



xa = £^e* (63) 

=0 V s 



1= 
n+l 



i=0 v« 



Using eq. (|63"D , the Mandelstam invariants can be written as [^3 



n+l 

s = x a xbs = E hikj\_e yi ~ Vi 

i,j=0 

n+l 

J ai = -2 Pa • fct = - ^ ^iie-^ (64) 

n+l 

s bi = -2p 6 • fci = - X ki±k jX e Vi ~ yj 

j=0 



Sij 



2ki ■ kj = 2k il _k j ^ L [cosh(y f - y 3 ) - cos(fa - fa)] . 



Eq.(|64D generalizes eq.flH]). As in that case, we note that the Mandelstam invariants 
depend only on rapidity differences, which again reflects the boost invariance of the 
center-of-mass frame with respect to the lab frame. 

Next, we assume to work in the kinematic region where the outgoing partons are 
strongly ordered in rapidity and have comparable transverse momentum, of size k±, 

Vo > yi > ••• > Vn+i] kn_ ~ kx, (65) 
From eq.(|64[) and (|65|), we obtain 

s ^ fto+ftn+iie" +1 , 
27 



s m ^ -h±h±e vo v \ 



(66) 



then eq.(^) can also be written as [|4j], p5| 



s > % > ki ± , 



(67) 



n 

2=0 



i=i 



The kinematic region where eq. (p5|) , or eq.(^), is valid is called the multiregge kinematics 
22|| , [Z3j|-||27|. The parton momentum fractions (^3]) assume here the simple form 



x A — 



B 



e 

n+l± 



k 



(68) 



which trivially generalizes eq . (|55D . Then we introduce the vector qi = p a — /co ; which 
labels the first momentum exchanged in the t channel (Fig. |6|). Using eq.(^), (|l^) and 
(B3), and keeping only the leading terms, 



q 1 ~ (k 1± e y \ -k ±e~ 



(69) 



Squaring qi and retaining only the leading term, 



q\ ~ 



(70) 



i.e. only the transverse degrees of freedom are relevant in the momentum transfer q\. 
The analysis can be repeated for the second momentum exchanged in the i channel, 
<?2 = Qi — ki, with the same conclusions, and so on for the following vectors qi, i.e. in 
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multiregge kinematics only the transverse components are relevant to parametrize the 
propagators of the gluons exchanged in the t channel, 

ii = q- -d±- (71) 

3.2 Three-parton production in multiregge kinematics 

First, we compute the simplest amplitude in the multiregge kinematics, the 2 — > 3-gluons 
amplitude (Fig. |). As discussed in sect. |2]^, we may consider just diagrams with gluon 
exchange in the t channel, as long as we work in a physical gauge. We consider first 
the diagram of Fig. |]a. Using the decomposition (^) for the gluon propagators in the t 
channel, the corresponding amplitude is 




Figure 4: 2 —>■ 3-gluons Feynman diagrams, with gluon k\ emitted a) from the gluon 
exchanged in the t channel; b) in the initial-state bremsstrahlung from the upper line; 
c) in the final-state bremsstrahlung from the upper line. We have neglected the helicity 
labels, and we have omitted the diagrams with bremsstrahlung emission from the lower 
line. 
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• g s f clC2dl [g vlV2 {qi + q2T 1 +g» 1V2 {-q2 + Piy i -g v ^ 1 {pi + <iiT 2 } 
■ (-*)P u b 2 ^(gsf bd2C2 g, b » 2 ), (72) 

st 2 v 7 

where we have implied physical polarizations for the external gluons. Then, we introduce 
the vector C^(qi,q2), which describes the contractions on the g^&i-ghions vertex, 

C"* (<&,<&) = - §P 7 [g UlU2 (qi + q2r+g" lU2 (-q2+Pir-g UllM1 (pi + qi) U2 }p u b 2 

* (Ql + Q2)T- § fp" b 1 + § fK 1 , (73) 

where we have used the approximate form (|66| ) of the Mandelstam invariants, and the 
decomposition fl39|) of (qi + 52)^, 

(qi + <&)" * (qi + Q2)l + ^rf - ^tf , (74) 

s s 

with = (0,0, qj_), such that q± ■ q± = —q±- Replacing eq. fl73]) into eq.([72|), the 
amplitude for the emission of gluon k\ along the {-channel gluon is 

■ Ug s r d ^C^( gi ,q 2 )) I (75) 

■ (wsf bd2C2 g, b , 2 ). 

The calculation of the amplitudes for the initial and final bremsstrahlung emissions 
(Fig. £|b,c) from the upper line yields respectively 

^ M ( ■i\ abd ° dld2 ~ °S n 3 jqcidi JCid c 2 fbd 2 c 2 fa 

uu r^ a /i t /iowfi 2 — zt, g s j j j ~ ? y v-a y M m 1 \ iyj ) 

i M( f\ ab (k)did2 95 3 fcidodi faeiC2 fftd 2 c 2 _^a_ 

Sal £2 
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Reordering the color with the Jacobi identities 



jac\ci jdad\c\ _|_ ydociC2 jd,\ac\ _|_ jd\c\C2 j 



cfociC2 fd\ac\ 



0. 



(77) 



the sum of the amplitudes for the bremsstrahlung emission from the upper line gives 



z(M(i) + M(/)) 



dz 



2is (ig s f 



ad{) c\ 



9/J.a 



/-'() 



u 



V s al I t% 



(78) 



(ig s f 



bd 2 c 2 



9n 



Adding the amplitudes for the emission of the gluon from the t channel fl75|) , the 
bremsstrahlung emission from the upper line ([78]), and the bremsstrahlung emission 
from the lower line (whose calculation is similar to the one of eq. (|76|) ), we obtain the 
2 — > 3-gluons amplitude in the multiregge kinematics, 

1 



j]\/f abd d l d 2 



where 



(?i + q 2 T± 



2is (ig. f aa ^ g, a 

(ig s f cldlC2 C% 
(W M2C2 W 2 ) 



Sal „ h \ Ml 

— + 2— pb 

s s bx 



1,Q2)) j- 



(79) 



Sbl 



+ 2^ 



-Sal 



Pa 



(80) 



is the non-local effective Lipatov vertex JMJ , which summarizes the insertion of the third 
gluon along the t channel, and as a bremsstrahlung gluon. The Lipatov vertex is gauge 
invariant, indeed replacing the physical polarization with the longitudinal one, the Ward 
identity 

c**Gzi,g2)(fci)„ = o, (8i) 
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holds. Now we want to show that the amplitude ([73) yields indeed ln(s/|t|) corrections 
to the Born- level two-parton production considered in sect. |2]4]. We square the amplitude 
(ff9|) and sum over helicities and colors. The sum over the helicity of gluon k\ may be 
performed using eq.(E3) since the Lipatov vertex is gauge invariant (pi). 



C»(qi,q 2 )-C,( qi ,q 2 ) = 4%^, (82) 

where we have used the kinematic constraint s a i<s&i = k 2 ± s from eq.(^). The sum over 
the helicities of gluons k and ki is performed using eq.([4^). By using then eq.([7j]) and 
([TT|) the square of amplitude (fT9|), summed (averaged) over final (initial) helicities and 
colors, is 

2 lQN?g 6 s s 2 

1^99^9991 - Tr2 7 p p p • 

Eq.(|33"D may be also obtained by taking the exact (square) matrix elements for three- 
gluon production in gluon-gluon scattering HlPf , 



\M„ g \ 2 = 4 (vra s iV c ) 3 J] 4 £ 1 , (84) 

i>j [a,0,l,2,6]' S aOS lSl2S2bSab 

with i,j = a,0,l,2,b, and with the second sum over the noncyclic permutations of the set 
[a,0,l,2,b]. In the strong rapidity ordering (|65|) and (|66|), eq.(Q) reduces to eq. (|83|) [ |45]j . 



Next, we must examine how the phase space for three-parton production, 

= f[ / (2vr) 4 *V fci) (85) 

transforms in multiregge kinematics. As we have seen in eq.(BBI), light-cone momentum 



conservation at large y (|55| ) simply generalizes to the multiregge kinematics. Accordingly, 



eq . (|57D also applies to three partons, and using it in eq.(^) the phase space becomes 

r 1 d 2 k 0± (dy 1 d 2 k 1± \ d 2 k 2 ± . /^;2 ^1 ioa\ 
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8=0 



which straightforwardly generalizes eq. (|58|) . The cross section for three-gluon production 
is then 



gg^ggg 



— Vi\M I 2 
2£ ' 3 \ lvl gg->ggg\ > 



(87) 



and using eq.(|83|) and (|86[ ) and the transverse- momentum conservation it becomes 



da 



gg^ggg 



y 



(88) 



47T kl ± kl L (kl ± + + 2k 0± k 2 ± cos 0) ' 

where we have performed the integration of the rapidity of gluon k± over the range 
y = yo — 2/2, and where is the azimuthal angle between the transverse momenta ko± 
and &2±- Since y = ln^/Arf), eq.(|38"D shows the logarithmic enhancement of the three- 
parton production with respect to the two-parton production considered in sect.[2~l. 
It also shows that the intermediate gluon is produced with equal probability over the 
allowed range in rapidity, i.e. eq.(^) is insensitive to the position in rapidity of the 
intermediate gluon. This approximation has important phenomenological consequences, 



and will be modified by subleading corrections, as we will see in sect. [O . 
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Figure 5: One-loop corrections to the elastic amplitude gg — > gg, with t-channel exchange 
of two gluons a) in the s-channel physical region and b) in the -u-channel physical region. 
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3.3 Virtual radiative corrections 



Next, we consider the one-loop corrections to the subprocess gg — > gg. In the large-y 
limit the diagrams with t-channel exchange of two gluons in the s-channel physical region 
(Fig.[|a) and in the w-channel physical region (Fig.[5]b) contribute. Let us consider the 
amplitude of Fig.[|a and parametrize the momentum transfer q 2 = i by the vector q ~ 
(0, 0, q_i_), where following the discussion at the end of sect. 3.1 only the transverse degrees 



of freedom are considered. Using the decomposition (39) for the gluon propagators in 



the t channel, the leading contribution to the amplitude of Fig.||a is 

iM ?Z^« - 2I 3 gt r "c r "a>o f w< cfb > c i gwg ^ bi I} (89) 

with / the integral over the gluon propagators, 

dad/3d 2 k± 1 1 



(27r) 4 a(3s - kj_ + ie -(1 - a)/3s - kj_ + ie 
1 1 



a(l + (3)s — kj_ + ie a(3s — (q — k)\ + ie ' 
where we have decomposed the momentum k on the light cone as 



(90) 



F = apg + prf + kj_, with < a, (3 < 1 , (91) 

and d 4 k = (s/2)dad/3d 2 k± the phase-space measure. Eq.(^) is called Sudakov decom- 
position. In eq.(|9~0"D all the propagators, but the second, have a pole in the lower complex 
half-plane of /3, so the integral over (3 may be easily performed in the upper half-plane, 
and we obtain 



/ 



dotd 2 k\ 
(2tt) 3 as - k\ k\ (q - kf 
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The integral of a is then logarithmic over the range {k\/s) < a < 1. Performing it and 
substituting eq. (|92|) back into eq.(^), the amplitude of Fig.|]a becomes 

where the adimensional function a(t) collects the loop transverse-momentum integra- 
tions, 

<*) J^ijtobjr (94) 

The integral in a(t) may be evaluated by introducing an infrared cutoff \x 

a s N c q\ 
In — 7, 

which shows that the amplitude ( |9~3"D is doubly logarithmic divergent. The amplitude of 
Fig.|||b is obtained by crossing the channels s and u in eq . (|93|) . Using then u = —s — i, 
and summing the diagrams of Fig|| the amplitude for the exchange of two gluons in the 
t channel is, to leading order in s/t, 



a(t>-^ln^, (95) 



iM z'z^ b . - ^! wwa(i) (96) 



The color dependence may be decomposed as a sum over the SU (3) representations in the 
Kronecker product of the adjoint representations, (8® 8), for the two gluons exchanged 
in the t channel, 

iM tZ a „y = ig^sT* Y,Pw(T)A T (s,t), (97) 

T 

where we have also singled out the helicity dependence. A T (s,t) are then scalar ampli- 
tudes, and ' (T) are color projectors 

P${T) P b J(T') = p£{T)6tt> • (98) 
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For the exchange in the t channel of a representation in the antisymmetric part of the 
product, (8 © 8)^ = 8 A © 10 © 10, or of a singlet, which is in the symmetric part of the 



product, (8 © 8)s = 1 © 8 S © 27, the color projectors are [06 , [47 



1 

1 



Pw (I) - N 2 _ i ^ ^ w ' 

aca' rbcb' 



P b a b H8 A ) = —r a f bc \ (99) 



1 



Nr. 

1 

c 



P b V(m®W) = -(5 ab 5 a , bl -5 abl 5 a , b )-—r a f 



The projectors (|9~9"D are normalized in order to satisfy the projection rule (|9~8"D . But for 
the singlet which is going to be of later use, we do not report explicitly the projectors 
of (8 © 8)5, since the leading contribution of it to eq. (|96|) cancels out because of the 
projectors parity under s <-> u crossing, 

P b r b '(T) = (-1) T P$'(T), (100) 

with 

(-1) T = 1 - 1 for ^®8U 

1 l) \+l for (8©8) s , 1 j 

The contraction of the structure constants in eq. ( |96| ) with P h a b °;' (10 © 10) from eq. (|99|) 
cancels out exactly, so the only contribution to leading s/t to eq.(|96|) comes from the octet 
in (8 ©8) a- Computing then the contraction of the structure constants with Pfffi (8 A ), 



eq. (|9~6"D becomes, 



iM Z'Z^ * 1 f ada 'f Mb ' ln Z| «(*) » ( 102 ) 

which yields the 0(a a ) virtual correction to eq. 
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We note that between the subleading one-loop corrections we have neglected there 
are the self-energy and vertex-correction diagrams which determine the running of the 
coupling constant. Accordingly, in a leading logarithmic treatment of the radiative 
corrections a s must be regarded as fixed. 

3.4 Multiparton production in multiregge kinematics 



In this section, we outline how the results of sect .[372] and |3.3| are extended to higher 
orders. In multiregge kinematics, eq. (f79D generalizes to the emission of an arbitrary 
number of gluons f25fl , namely the tree- level multigluon amplitude preserves the ladder 
structure (Fig. |a), and for the production of n + 2 gluons it is given by 
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Figure 6: Multigluon amplitude in multiregge kinematics a) at tree level and b) with the 
virtual radiative corrections, represented by the thicker gluon line. The blobs remind 
that non-local effective Lipatov vertices are used for the gluon emissions along the ladder. 
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(ig s r d ^c^( qi ,q 2 )) i 



(103) 



(ig s r d ^C^{q n) q n+l )) J- 

"n-\ 

(ig s / M "+ lC ™+ 1 9ilh 



where the Lipatov vertex (p0|) is 



(ft + 



(104) 



A proof of the validity of eg. (|103|) may be found in sect. 4 of ref.[[H 



Following ref. ||24|| , we make the ansatz that the leading logarithmic approximation 
(LLA) in ln(s/t) to the virtual radiative corrections, to all orders in a s , is obtained by 
replacing the propagator for the i th gluon with 



a(U) 



—e 
U 



(105) 



with a(ti) as in eq. fl94|) . Inserting eq. (|105|) into eq.(|40|) we obtain the amplitude for the 
subprocess gg — > gg, with the leading contribution of the virtual corrections to all orders 
in a s , 



iM, 



aba'b' 



-8ma s r' c g, atlal je a ^^^ f bb ' c 



9n, 



b Mb' • 



(106) 



In eq. ( |102| ) we have reproduced the calculation of the 0(a 2 s ) term of eq. (|106| ) , the 0(a? s ) 
term has been computed in ref. |24|, and in ref.[[25|] eq.( |106| ) has been conjectured to be 
valid to all orders in a s . This entails that the leading contribution of the virtual radiative 
corrections to any order in a s has the color-octet structure of one-gluon exchange. Thus, 
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in the decomposition (|97p of eq.( |106| ) the only scalar amplitude that contributes to 
leading s/t is, 




Replacing eq. (|95|) into eq. ( |105| ) we note that the argument of the exponential factor on 
the right-hand side is a negative product of a collinear logarithm and a logarithm of 
type ln(s/t). This is an example of Sudakov form factor. Accordingly the amplitude 
( |106|) and the related production rate are infrared sensitive and vanish as we take fi — > 0. 
This infrared behavior is general in gauge theories and was first noticed in QED [ ff8|j , 
where the probability for having a scattering process without the emission, or with the 
emission of a finite number, of soft photons vanishes. It is only after the inclusion of the 
real radiative corrections that a cross section is well defined and finite order by order in 
perturbation theory 

Inserting then eq. (|105|) into eq. (|103| ), we obtain the multigluon amplitude for the pro- 
duction of n + 2 gluons, with the leading contribution of the virtual radiative corrections 
to all orders in a, 



j J\fCibd ...d n+ i 



) 



g a (*l) (2/0-2/1) 



tl 




(108) 




(ig s / M «+i c "+i 9ltb 



Eq. ( |108j ) has been verified for three-gluon production at one loop in ref. [p5|] . A proof 
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of its validity to all orders in a s may be found in the appendix of ref. |[46|| . We will 
follow, though, the original approach of ref. |25fl , namely we take eq. (|108| ) as a corollary 



of the ansatz ( |105| ) and we prove its self-consistency by using eq.( |108|) to derive the 
elastic amplitude with the virtual radiative corrections in LLA to all orders in a s , and 
by showing that it coincides with eq. (|106|) . 



3.5 Partial wawe amplitudes 

In this section we derive dispersion relations which will let us reconstruct the elastic 
amplitude with all the virtual radiative corrections ( |106|) . 



First, we decompose the scalar amplitude A T (s,t) introduced in eq.(|97|)in t-channel 
partial wawe amplitudes 

A T (s, t) = £ (21 + 1) Af(s, t) Pi{zt), (109) 
i 

where / is the angular momentum, Pi(z) are Legendre polynomials and zt = — cos 6t is 
the scattering angle in the t-channel physical region, whose Mandelstam invariants are 
obtained from eq. ( |I~9"D by crossing the channels s and t, 

8 = -^(1-COS^), (HO) 

u = --(l + cos^). 
Using the invariance of amplitude (^7|) under the crossing symmetry s *->■ u, 
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and the parity of the projectors ( |100| ), we obtain the parity of the scalar amplitudes 
under s <-» u crossing 

A T (—z t ,i) — (—l) T A T (zt,i) , (112) 

with the change of sign of the scattering angle, z t *-> —z t under s <-> u crossing, as seen 
from eq. (|110|) . Next, we write the elastic amplitude (p7f) through a dispersion relation, 
i.e. as an integral over the amplitude singularities, which are branch cuts over the real 
axis of the s complex plane, — i < s < oo for the physical s channel, and — i < u < oo 



for the physical u channel [50 



At**- f° ds ' Di sc^(s',t) r ds ' DiscA(s',t) 

A[8,t) ~ LootoU s'-s + Li 2ri s'-s ' [1U} 

with 

Disc A(s', i) = A(s' + ie, i) - A(s' - ie, i) , (114) 
and where we have used u = — s — t, and i < 0. Since, from eq. (|110|) , 



zt=~u+jy (us) 

the dispersion relation (|113|) may be written as an integral over the complex plane of z t , 

WA ~ x r- 1 dz' f Disc A(z't,t) r°° dzL Disc A(z' f , t) 

A(s,t)= — , w ' ; + / -4 , V * ■ 116 

j-oo 2m z t — z t Ji 2m z t — z t 



We note that in the i-channel physical region (|110|) , where z t is the scattering angle, 
— 1 < Zt < 1, while the singularities are at Zt < — 1 and Zt > 1, where the t channel is 
unphysical ( |116| ). So at — 1 < z t < 1 we may invert the partial- wawe expansion ( |109| ), 
to obtain the amplitude for the I th wawe 

Af{s, t) = l - £ dzt Pt(zt) A T (s, t) . (117) 
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Introducing the Legendre function, associated to the Legendre polynomial 

l r 1 dz 

2 



Qi(z') 



1 z — z 



-Pi(z), 



(118) 



replacing the dispersion relation (|116|) into eq.( 117 ), and using the parity under s <->• u 
crossing, 



Qi(-zt) 
Disc A T (-z u i) 



the amplitude for the I wawe becomes 



Aj(s,t) = l + (-l) 



\l+T 



-i) m Qi(zt) , 

-1) T+1 DiscA T (^,t) 



00 d 7 1 

— Qt(z') Disc A T (z',t) 

1 Z7T? 



(119) 



(120) 



Next, we introduce the Sommerfeld- Watson representation of the amplitude A T (s,i) in 



the complex plane of the angular momentum / |5T 



A T (s,i) = f.(2l + l)Af(s,t) 

J 0—ioo Zt 



sin7rZ 



(121) 



where the integral is done along a path parallel to the imaginary axis, to the right of 
all the singularities. All that we have considered so far in this section describes general 
properties of the analiticity of the scattering amplitude. Now, we consider the high- 
energy limit s —>■ 00 at fixed t. From eq. (|115|) we have, 

(122) 

and we take the asymptotic value for large z t of the Legendre polynomials and their 



2s 

Z t ^-j 



associated functions 52 



Qi(z) 



1 rf/ + i 



r(/ + i) 
rq + i) 
r(/ + §) 

42 



(2z)-^ 



(123) 



Replacing the amplitude for the I th wawe ( |120| ) into the Sommerfeld- Watson representa- 
tion of the amplitude (|121|) , and using the asymptotics ( |123| ), we obtain 

A T (s, t) = -i- f 5+lQ ° dl ( " 1)r e^(t) , (124) 

47T J8~ioo Sin 7Tt 

where Tf{t) is the Laplace transform of the discontinuity of the amplitude 

= / rf?/e^DiscA T (z t ,t), (125) 
Jo 

with y = In (z t /2). 

3.6 The BFKL equation 

Our goal is now to use the s-channel unitarity and the multigluon amplitude with ex- 



change of one gluon in the t channel (|108|) to evaluate the discontinuity of the amplitude 
for exchange of two gluons in the t channel in a physical gauge. From that by using the 
dispersion relations computed above the elastic amplitude with all the virtual radiative 
corrections ( |106j ) will be derived. 

In order to evaluate the discontinuity of the amplitude, we compute and sum over 
n the cut diagram of Fig.|7|. In evaluating the diagram at n + 1 loops, we must put the 
internal legs fcj on mass shell, by replacing the n + 2 internal propagators with the cut 
ones 



153 



-^2n5(k*). (126) 

Multiplying then the mass-shell factors (|126|) by the n + 1 loop-momentum integrals, we 
obtain the phase space for the production of n + 2 partons 

"•^S/^S w(t]l = .3/ w W5 ' h+B "S y (l27) 
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where we have inserted the identity / d A k n+ i5 4 {p a +Pb~ ^7=0 = 1; an d use d eq.(|14|). 
Now we need to consider the phase space in the multiregge kinematics. This has already 
been computed for three-parton production in eq. fl86|) . Since light-cone momentum con- 
servation in the multiregge kinematics (|68|) has the same form for three or more partons, 
eq. (|86D trivially generalizes to n + 2 partons, and the phase space ( [1271 ) becomes, 

1 d 2 kn / n A - A2r 



mii / tt dyid ki±\ d k n+1 ± /ri _ N2 r 2 / 1.2 



2s (2 



7T 



n 



\ 4tt (2 



7T 



2 7. /n+1 

(2vr) 2 5 2 E^a 

\i=0 / 



(27T) 



Then using eq.( |108| ) the discontinuity of the amplitude of Fig.|7| is 



(128) 



Pa a 



kn dn 



Pbb 



1 qi ci 
7 k a di 




< k 2 d 2 




< k n d n 




| Qn+l Cn+l 

) k n +i d n+1 





q - q2 c 2 



q - qn+i c n+1 
Pb-b' 



Figure 7: 0(a™ +2 ) contribution to the discontinuity of the elastic amplitude with two- 
gluon ladder exchange in the t channel. 
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(ig s f CndnCn+1 ) (ig s f <dnC '^) C^(q n ,q n+1 )(-g^JC^(q- q n , 



Q - Qn+l) 



3 a(f„ + i) (y n -y„+i) 



3 a(*n+l) (Vn-Vn+l) 



■ Ug s f bdn+1Cn+1 ) (ig s f c ' n+ldn+lb ' 

where t = q 2 is the momentum transfer in the elastic scattering and t[ 



- Qi) 



In eq.( |129| ) we have used the phase space (|128|) , and we have done the sum over the 
helicities of the gluons emitted from the helicity-conserving vertices using eq. ([42"D , and 
the ones over the n intermediate gluons using eq.([RJ) since the Lipatov vertices are gauge 
invariant flSlf). The Lipatov vertex C^(q — qi, q — qi+i) in eq. (|129|) may be obtained either 
by direct construction as in sect. |3^ or by simply inverting the sign of the momentum 
ki in eq.fllPID, 



C^(q-qi,q-q i+ x) 



V s s bi J 



Sb 



t' 



(130) 



Then using eq. ([7T|) , ( |104|) , (|130|) and the kinematic constraint s ai Sbi = kf ± s from eq.(|6 
with kf ± = (qi — qi+i) 2 j_, the contraction of the Lipatov vertices becomes 

\2 J2 



C^(q h q i+1 ) C^q-q^q- q i+l ) 



2 (g ~ Q i)±Qi+i± + {Q- Qi+i)±Qh 



(ft - Qi+i)± 



= -2fC(qi,q i+1 ) 



(131) 
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Next, using eq.([J7]) we consider the scalar part DiscA T (s, t) of the discontinuity ( pL 29| ) . 
The contraction of the color projectors (f^) with the structure constants in eq.( |129| ) 
yields the color factor C^ +2 , with 

r _{N C for the singlet, . . 

Ut ~ \ N c /2 for the octet, [i6Z) 

where we have used the Jacobi identity ( |77|) in the contraction with the octet projector. 
So the scalar part of the discontinuity (|129|) may be written using eq. (|131|) and ( |132| ) as, 



Disc^ (s ,t) = E(-^) n+2 /n$n 



\ (2tt) 2 



n=0 i=l j 

ri+1 i n 

■ 2 ^ II 717 e^-y^ +a ^ J] 2/C(g m , g m+1 ) , (133) 

i=l HH m=l 

where we have used the conservation of the transverse momentum and we have changed 
integration variables from the transverse components of the produced gluons to the ones 
of the gluons exchanged in the t channel. In eq. (|133|) we have n integrations over the 
rapidities of the gluons emitted within the ladder, while in the integrand there are n + 1 
rapidity differences. To disentangle the integrations, we take the Laplace transform 
( |125|) of eq. (|133|) with respect to the rapidity difference y = yo — y n +i = In (— s/i), and 



change the integrations over the rapidities y and yi, with % = 1, n to the ones over the 
rapidity differences — yi, with i — 1, ...,n + 1. Then we perform the integrations, 
with boundaries given by the strong rapidity ordering (j65|) , and we obtain 

oo . n+1 J2 _ 

= -2.t(4™ s ) 2 ciE/n^% 

n =cr j=l \ Z7T ) 

1 1 



t x t' x I -I- a(ti) - a(t[) 
(-2a s C T ) K,{q u q 2 ) 
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1 1 

W 2 l-l- a(i 2 ) ~ a(t' 2 ) 



(134) 



-2a s C T ) K{q n , q n+l ) 
1 1 



i n +li'n+l I ~ 1 — a (*n+l) — a \t'n+l) 

This may be written as a recursive relation, 

d 2 qi± 1 



T?{t) = -2it{^a s ) 2 C 2 T J 



//(ft,*), 



(135) 



where the function /; T (g, i) satisfies the integral equation, 



/ - 1 - a(t x ) - a(t[) f?( qi , t) = 1 - 2a s C T 



d 2 g 2 _L K.(q 1 ,q 2 ) 



f?{q 2 ,t), (136) 



with = —qf± and £■ = —(q — g^) 2 . Eq.( |136|) is the BFKL integral equation 
describing the gluon-ladder evolution in the LLA of ln(s/t). The function K(qi,q 2 ) 



defined as the contraction of the Lipatov vertices ( |131| ), describes the contribution of the 
real radiative corrections and forms the kernel of the BFKL equation. The contribution 
of the virtual radiative corrections a(i) appears on the left-hand side of eq.( |136| ). 

Next, we are going to solve the BFKL equation for the color-octet exchange. Replac- 
ing then eq. ( |132|) and the explicit form ( |94"D of a(i) into eq. (|136|) , this reduces to 

d 2 ki 1 



(l-l)fr\qi,t) = l-a s N c q{ 
which admits the solution 

fi oct (qi,i) = fT*(k,t) 



(2tt) 2 k\(q-kf_ 



(137) 



1 



I -I- a(t) 



(138) 
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The solution is unique, since eg. ( |136| ) is inhomogeneous, and it is constant with respect 
to the functional dependence from its first argument. Replacing eq.( |138| ) into the Laplace 
transform (|135| ) and using eq.(|94]) and (|132|) , we obtain 



rr\t) = -8vr 2 * N c a s a(t) . (139) 

I — 1 — a(t) 



Replacing eq. ( |139| ) into the Sommerfeld- Watson representation of amplitude ( |124 ) and 



using the parity factor ( |101| ), the scalar part of the amplitude for octet exchange becomes 

l+a(t) 



A° c \s, t) = 4vr N c a s -^L (l + e-©) [A ) . (140) 



where the integral over the complex plane of / has yielded a Regge pole at I = 1 + a(t), 
which then gives the Regge trajectory Since a(t = 0) = 0, the intercept is at 



1 = 1, which corresponds to a reggeized gluon [25], |51| 



Even though we started with a discontinuity (|129|) at 0(a 2 s ), that computes the 



radiative corrections due the insertion of a two-gluon ladder, the octet solution ( |140| ) is 
0{a s )i because the octet appears already at the one-gluon exchange level. 

Finally, we note that in the multiregge region (|65|) , where the collinear logarithms 
are not supposed to be large, a s In (gi//i 2 ) 1. Then from eq. (|95|) we obtain that 
7r|a(t)| 1. Taking this limit in eq. ( |140| ) , we find it in agreement with eq.( |lQ7| ), proving 
thereby the self-consistency of the ansatz ( |105| ) for eq.(p 

3.7 The pomeron solution 



We are now interested to study the singlet solution of the BFKL equation ( |136| ), since it 
can be related via the optical theorem to the total cross section for one-gluon exchange 
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in the t channel. Let us write first the function fi ing (qi,t) in differential form, 

d 2 k< 



(27T) 



;fi(qi,k,t). 



(141) 



Replacing it into eq. (|136|) and using eq. (|132|) , we obtain the BFKL equation for color- 
singlet exchange, 



/ - 1 - a(ti) - a(i[) fi(q u k, i) = 

d 2 q 2 j_ JC(q 1 ,q 2 ) 



(2ir) 2 6 2 ( qi -k) - 2a s N c J 



fi(q2,k,t) . 



(142) 



(2tt) 2 q 2 ± (q~q2)i 

Since by pomeron exchange it is meant conventionally the exchange of something with 
quantum numbers of the vacuum, the exchange of a two-gluon ladder in a color-singlet 
configuration defines the exchange of a perturbative QCD pomeron. 



The kernel K,(qi,q2) ( |131|) is regular as q%± — > oo, making the right-hand side of 
eq.( |142| ) ultraviolet finite. There are ultraviolet divergences as q 2 ± — > oo in the virtual 



radiative corrections (|95D in the left-hand side of eq. (|142|) , but they cancel in the singlet 
solution (|135| ). As for the infrared behavior of eq. ( |142j ), the kernel vanishes at q%± = 0, 
or at (q — q^)! = 0, where accordingly the right-hand side of eq.( 142 ) is regular. The 
kernel is singular at q\ = q 2 , however making the substitution 



d 2 q 2 



d 2 q 2 



fa 



(143) 



in the virtual radiative-corrections terms of eq. ( |142| ), it is easy to see that the singularity 
at q\ = q\ cancels between the virtual and real radiative corrections. Still, eq.( |142j ) has 
infrared divergences as q\ — > in the virtual radiative corrections. Besides, at qf ± = or 
(q — <?i)i — the kernel vanishes, and after using eq. ( |141| ), eq. ( |142j ) reduces to eq.( |137| ), 
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thus the two-gluon ladder reduces to a one-gluon ladder, which is infrared sensitive as 
we know from sect. 3. 4 and 3.6. However, it is possible to show that in particular cases, 
like for the exchange of a singlet gluon ladder between colorless objects, the solution of 
eq. (|142|) has no infrared divergences at all [54 . 



Since we are interested in the forward scattering amplitude and in the total cross 
section, we look for the solution of eq. ( |142| ) at q = 0. From eq.(^) and ( |131| ) we obtain 



£(<?l,92)|g=0 

Doing then the substitution of variables, 



(<?i - ?2)± 



2 - 



1 k 2 
57T q 1± 

the BFKL equation for singlet exchange at i = becomes 



(144) 



(145) 



(l-l)fl(qi,k) 



a s N c 

7T 2 



/ d 2 q 2 ± t 



(146) 



fi{q2,k)-$^fi{ qi ,k) 



2q 



2X 



This is a inhomogeneous integral equation with a self-adjoint kernel, given by the real 
radiative corrections to the discontinuity of the amplitude. As we have done for eq.( |142| ), 
in order to regulate the divergence at q\ = q\ we make the substitution ( |143| ) in the 
virtual radiative-corrections term on the right-hand side of eq.( 146|) (cf. Appendix A), 
then the BFKL equation at t — reads 



(I - 1) f t ( qi , k) = 



[ d 2 q 2 ± t - 



fi(Q2,k) 



Q2± + (Qi ~ 



(147) 

"2 /i(?l>*0 
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The homogeneous equation associated to eq. ( |147| ), may be written as 



(I- l)fi(qx,k) 
®sN c f l2 f fi(q 2 ,k) 



(148) 



7T Z 



fr2 J M 



i± 



(Qi - 92)1 ql± + (qi ~ 42): 



fi{qi, k) 



and admits a solution Fourier series 

fi{qi,k) = £ / cW^)e MAl " A) e m( 



6) , (149) 

with Ai = ln(g 2 /m 2 ), A = \n(k 2 /m 2 ), m 2 a scale factor, and <p\ — <fi the azimuthal angle 
between the vectors k± and q±±. Using then the integral representation of the 5-function, 



the inhomogeneous term in eq. ( |147| ) can be expanded as 



S 2 (qi - k) 



1 00 



dve iu{Xl - x) e in{01 



(150) 



(k 2 ± q 2 ± y/ 2 2tt 2 

Replacing eq.( |150| ) into eq.( |147| ) and calling u(u,n) the eigenvalue of the homogeneous 
equation (|148|) , we obtain the condition 

1 1 



(/ — 1) a(z/, n) 



(k 2 ± q 2 ± y/ 2 (2tt) 2 
which fixes the coefficient a(z/, n) in the solution ( |149| ) 



+ cj(z/, n)a(is, n) 



(151) 



fi(qi,k) 



l l 00 roo I 

hi) 2 (k 2 q 2 Y' 2 J^J-^ l-l-u[u.,n 



e iv(Ai-A) e «n(</>i- 



(152) 



(2n) 2 (klq 2 ± y 

To find the spectrum of eigenvalues, we replace the solution ( |152| ) into the homogeneous 



equation (|148|) , and obtain 

a s N c 



ui(u } n) 



71 



X 2 



2 Re / dx ■ 

1 — x 



-1 l rl 

— I dx h / dx 



2 / dx 

JO 1 — X 

+ f dx 



x Jo xa/1 + 4x 2 J Vx 2 + 4 



(153) 
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with 

The last three terms in eq.i 
of the T function 



x 



ql/ql for q\ < q\ 



(154) 



~ \ Qi/ql for q\ > q\ . 
cancel out, and introducing the logarithmic derivative 



d\nT(y) 
dy 



1 , a*" 1 - 1 

ax 

o x — 1 



7' 



(155) 



with 7 = 
becomes 



0.577215... the Euler-Mascheroni constant, the eigenvalue ( |153| ) 

.a„N K 



uj[y^ n) 



-2^^ Re 

7T 



1> 



n 



+ iu - V(l) 



(156) 



uj(y,n) is regular at both the integration limits (155), entailing that eq. ( |147| ) is regular 
at q\ = 0, q\ = qf, and q\ — > oo, i.e. that all the infrared and ultraviolet divergences 
cancel out. 



Next, we use the solution ( |152j ) in the Laplace transform of the discontinuity of the 
amplitude ( |135| ). The right-hand side of eq.( |135| ) is proportional to t, so it seems to 
vanish. However, the factor i was due to replacing s = —ie yo ~ y " +1 into eq.( |133|) , and for 
the t = case we should have rather done the substitution s = kj_e yo ~ Vn+1 , with k\ as 
in eq. (|B3D . Doing so, and using eq.( |132| ), ( fL41| ) and ( |145| ), the Laplace transform of the 
discontinuity for singlet exchange becomes 

1 



J^(t = 0) = mk\a 2 s N 2 c f(fka±cPk b±w -^-f l {k a ,k b ) 



(157) 



with fi(k a ,kb) as in eq.( |152| ), with eigenvalue ( P-56| ). Before taking the inverse Laplace 
transform of eq.( 157 ). 



Disc A(s,t) 



dl_ 

2vri 



158) 
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it is convenient to examine the eigenvalue (|156|) to determine where the leading singu- 
larity at large y is. The leading contribution to n) comes from n = v = (cf. 
Appendix B), and for small v the eigenvalue admits the expansion (|189| ) 



uj(u,n = 0) = (21n2-7C(3)z/ 2 + ...) . (159) 

7j- V / 

Replacing it into eq.( |152j ), and evaluating the integral over the complex plane of z/, we 
obtain the leading contribution to the solution for the singlet, 

f,(k kh) ~ — - - -^o| ln(fc2 /fcg ) 

b) ~ W (ki x ki ± y/* m - 1 - a)} 1/2 ■ [ } 



with 



A = 41n2«^ ; B =14C(3)^; (l^_A) U \ (161) 

7T 7T V D I 



Replacing it into the discontinuity ( |158| ), we see that in the complex plane of I the leading 
singularity of the pomeron solution is a branch cut, extending from — oo up to / = 1 + A 



Now, let us go back to the full solution ( |152j) . Taking the inverse Laplace trans- 
form ( |158| ) of eq.( |157| ) and using eq.( |152[ ), the discontinuity of the amplitude for singlet 



exchange is 



DiscA s ^(s, t = 0) = 16* s a] N 2 f d 2 q a± d 2 q bl _ -^f(q a , q b , y) , (162) 



with y = hi(s/kx) and, 



/(*,*,*) = T^wh^ £ e^jy^e^J^, ( 163 ) 



(2vr) 2 (q 2 a± q 2 b± )^ n 
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with the azimuthal angle between the vectors q a ± and q b j_, corresponding to the first 
and the last gluon momenta exchanged on the ladder (Fig||). f(q a ,qb,y) may be also 
written as the inverse Laplace transform of the solution ( |152|) with respect to ou — I — 1. 

Next, we consider unitarity and the optical theorem to relate the total cross section 
for gluon-gluon scattering with exchange of a one-gluon ladder (Fig.^D) to the forward 
amplitude for gluon-gluon elastic scattering with exchange of a two-gluon ladder (Fig .[7]), 

1 r , + 1 



°tot 



— f dVM M 1 " = Disc \iM(s, i = 0)1 , (164) 

2s J 2s 1 J 



with V the phase space of eq. ( |127| ), and where we have used eq. ( |114|) . M(s,i = 0) is 
the forward amplitude summed (averaged) over final (initial) colors and helicities. From 
eq. (|97D it can be written as 

Using then eq.( |162| ) and (|165|) , the total cross section for gluon-gluon scattering ( |164j) in 
multiregge kinematics is given by 

crtot = AT2 c : a 2 s / d 2 k al _ d 2 k b ± 2 2 f(k a , h, y) , (166) 

with k a± = -q a± and k bl _ = q b± , as in Fig.g f(k a ,k b ,y) is given by eq. ([T63|) with 
= + ii the azimuthal angle between the vectors k a ±_ and k b ±, and y = In (s/k a j_k b j_) 
the rapidity interval between gluons k a and k b . 

As noted after eq. (|156|) , the solution ( |152| ) or (|163|) of the BFKL equation for singlet 
exchange at i = has no infrared and ultraviolet divergences. Namely, the doubly 
logarithmic enhancements of the virtual and real radiative corrections cancel out, leaving 

54 



Pa ka 


Pa' 


i 

o 

Pb o 




O 
O 

o 








qb 

k b 


o 
o 
o 
o 







Figure 8: Two-gluon production with exchange of a BFKL ladder, represented schemat- 
ically by the double-lined square. 

as a residue single logarithms. Accordingly, the expansion of eq. (|163|) is finite order by 
order in a s . This is due to the optical theorem ( |164| ) which relates the forward scattering 
amplitude to the total cross section, for which all the infrared and ultraviolet logarithmic 
divergences in the radiative corrections to the Born amplitude cancel out, order by order 
in a s . Besides, as discussed at the end of sect. 3.3 the coupling constant a s in eq. ( |1 47| ) 
must be regarded as fixed, thus the total cross section (|166| ) does not depend on a 
renormalization scale Finally, we note that because of the leading singularity ( |160| ) 
at I = 1 + A, the growth of the total cross section, otot = s A , violates the Froissart 
unitarity bound a to t oc In 2 s |51fl . 



''Attempts to introduce by hand a running coupling constant in eq.(147) have been made ]55| ], but a 
running a s makes cq.(147), and so the total cross section (166), infrared sensitive. 
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3.8 The total parton cross section 



Summarizing the conclusions of sect . 577 , the parton cross section for the production of 
two gluons, resummed to all orders of a s In (s/k\), is obtained from eq. (|163|) and ( P-66Q , 



da, 



99 



C A a s 



k 2 



f(k a , k b ,y) 



C A a s 



k 2 



(167) 



d 2 k a ±d 2 kb± 

with y = y a — y b and f(k a ,k b ,y) the inverse Laplace transform of the singlet solution 

/,(*«,**) o, 

f(k a ,h,y) 



du 



2m 



e" v fi(k a ,k b ) 



(168) 



with u) — I — 1. Eq. (|167 ) is written as a convolution of the singlet solution (|168|) with the 
parton-production vertices on each side of the rapidity interval. The structure of eq.( |167| ) 
is generic and valid for any scattering process in multiregge kinematics (cf. Appendix 
C). 



First, we consider the limit a s y — > in eq. (|167|) , for which all the real and virtual 



radiative corrections vanish, and the singlet solution (|168|) reduces to the inhomogeneous 
term of the BFKL equation, 



f{k a ,h,y)\o{a°) = ${kl± - kl ± ) 5(<f) - it) 



(169) 



i.e. at the Born level only two partons are produced, and they are balanced in k± and 
back-to-back in 0, and eq. (|167|) reduces to the Born parton cross section 



By integrating eq.( |167f ) over the azimuthal angle only the n = contribution to the 
eigenvalue (|156|) survives 



da, 



mi 



dk 2 al dk 2 ^ 



(170) 
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In the asymptotic limit of large rapidities, we may approximate the eigenvalue ( |156| ) 
with its expansion (|159|) at small z/, and perform the integral over v by the saddle-point 
method, 

da gg nC 2 A a 2 s / 1b 2 (A£l/*£l^ 



dP a± dP b± Akl L kl L y/nBy 6XP [ ABy ) ' 
with A and B given in eq. ( |161| ). Eq.( |171| ) shows a diffusion pattern, namely it has 



the form of a Gaussian distribution in hi(A;^ ± //^ ± ) with the peak positioned where the 
partons are balanced in k± and with a width growing with y. This is not accidental, 
since the BFKL equation may be written as a diffusion equation, with diffusion rate 
HklJH±) ~ V 1/2 (cf- Appendix D). 

Let us go back to eq. ( |170| ), and perform the integration over the parton transverse 
momenta above a cutoff p± m m- As noted at the end of sect. 3.3 and |3.7|, in the BFKL 



formalism the coupling constant is fixed. Accordingly, for all of the partons produced 
in the gluon ladder the coupling must be evaluated at a single scale. Typical scales are 
the transverse momenta k a ± and kb± of the partons with respect to which the ladder 
is being resummed, and their cutoff p_i_ m i n . Choosing as a scale any combination of 
them is theoretically equally valid, and anyway we expect from eq.( |171|) that different 
choices do not make a large difference within the BFKL approximation. Here we choose 
a s = a s (pi min ) plf , since it allows us to perform analitically the integrations over the 
transverse momenta, thus obtaining the total parton cross section, 

a « = tr L l dp ^r^- (172) 



Approximating then the eigenvalue (|156|) with its expansion ( 159|) at small u, and per 
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forming the integral over v, we obtain the large-?/ asymptotics of the total cross section, 

°gg = 7T^ — ~ / ( 173 ) 
2pl miD ,/7C(3)iV c a s y/2 



Eq. (|171j) and ( 173 ) show the exponential growth in y, corresponding to a power-like 



growth in s, characteristic of the pomeron solution ( |160| ). Dividing eq.( |172| ) or ( |173| ) 
by eq. (|54]) , with t integrated over the transverse- momentum cutoff Pj_ min , we obtain the 
.fT-factor, i.e. the growth rate of the total cross section due to the radiative corrections, 

1 poo pU>(i',n=0)y p 4ln2Nca s y/iT 

K=— dv— — ~ . (174) 

with y = \ia(s / p 2 ±min ) . In fig.|^ we plot the K-factor as a function of the parton center-of- 
mass energy \/I at different values of the transverse-momentum cutoff pj_ m i n . We scale 
a s = a s (pi m i n ) from a s (m 2 z ) = 0.12 using the one-loop evolution with five flavors. Even 
though the power of s in the K-factor is not very large (at p± m i n = 20 GeV, the power 
is ~ 0.4), the if-factor grows quickly [|. 

4 Two-jet production at large rapidities 



In this section inclusive jet production is considered as an experimental signature of 
the BFKL theory discussed in sect. 3. The original proposal of Mueller and Navelet 
[ 2j| of linking the growth rate of the inclusive two-jet production to the i^-factor of 



the total parton cross section is recalled, and its modifications to fixed-energy colliders 



l We note that the dashed curves approach the solid ones on the plot logarithmic scale, since the 
relative error tends to zero at large s, however on a linear-scale plot the distance between the curves 
increases, since the absolute error increases at large s [p3[. 
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Figure 9: K-factor as a function of at different values (in GeV) of p± m i n . The solid 
curves are obtained by doing numerically the integral over v in eq.( |174|) , while the dashed 
curves correspond to the asymptotic expression on the right-hand side of eq. (|174|) . 



are discussed; finally, the jet-jet decorrelation in transverse momentum is proposed as 
a probe of the BFKL dynamics, and an assessment is made of the phenomenological 
importance of higher-order and next-to-leading logarithmic corrections to the BFKL 
formalism. 
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4.1 The Mueller-Navelet jets 



In order to analyze two-jet production in a way that resembles the configuration assumed 
in the BFKL theory, we select all the jets in the event above a transverse momentum 
cutoff p± m i n , and rank them by their rapidity, i.e. we tag the two jets with the largest 
and smallest rapidity and observe the distributions as a function of these two tagging 



jets (in the standard hadronic-jet analysis jets are ranked by their transverse energy [11 
|l||])- Since the jet production is inclusive, the distributions are affected by the hadronic 
activity in the rapidity interval y between the tagging jets, whether or not these hadrons 
pass the jet-selection criteria. 

In order to obtain the inclusive two-jet production we must convolute eq.( |167| ) with 
parton densities, whose dependence on the parton momentum fractions is given by 
eq. (|68D , with ko = k a and k n+ \ = k^. Our goal is to examine the parton process ( |167| ) 
at growing values of y, to check if indeed it shows the exponential growth in y of the 
pomeron solution. At the same time it would be convenient to keep the x's fixed, in order 
to disentangle the eventual dynamical rise due to eq.( p.67] ) from kinematic variations in- 



duced by the parton densities. Thus, Mueller and Navelet have proposed to measure 
two-jet production at fixed x's and different values of y. This implies that the hadron 
center-of-mass energy yfs has to be risen proportionally with the parton center-of-mass 
energy \f§, i.e. with the exponential of the rapidity interval y between the tagging jets. 

We need a factorization formula which relates the parton cross section ( |167[ ) to jet 
production. The factorization formula for the two-jet production cross section in terms 
of the jet rapidities and transverse momenta has been given in general in eq.(pO[). The 
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parton cross section for the production of n + 2 partons is 

1 



dan = j-.dVn^lMijl 2 , (175) 



with the differential phase space (|127 ) for the production of n + 2 partons. Using the 



light-cone momentum conservation in the exact kinematics (|B"3D, and eq.( |175[ ), the two-jet 
production cross section (p0|) becomes 
da 



Of 2 fc^.L G? 2 fc fo _L C?2/ a G?2/ fo 

L y (11 y 47r(27r)2 j L /i ) x B fj/B{xB, /i ) j^tj 5 ( X, &u 



(176) 



n=0 ^ \/=l ,/ V"V / ij iu " ° \i=o 

where the n = term reproduces eq. (|22"D . In the large-y limit the phase space for the 
production of n + 2 partons is given by eq. (|128|) , where the approximate kinematics 



(|68|) is used, and the parton momentum fractions fix the jet rapidities (cf. eq.([57|)). 
Using eq.( |128| ) and ( |175| ) into eq.( |176| ) we obtain the two-jet production cross section in 



multiregge kinematics, 



</ e //(4^ 2 )4/e // (4^ 2 )^7%7-, (177) 



d 2 k a ±d 2 k b ±dy a dy b ' ' d 2 k a ±d 2 k b ± 

with the parton cross section as in eq.( |167| ). In eq.( |177| ) we have used again the effective 
parton density (|53|) since the leading contribution to two-jet production in multiregge 
kinematics always comes from the exchange of a gluon ladder in the i channel, the only 
difference between subprocesses with initial-state quarks or gluons being the different 
color strength in the jet-production vertices. From eq.(^) and (|177|) we can now easily 



derive the factorization formula at fixed parton momentum fractions, 

feff(A, A* 3 ) feff(A, f) M f» • (178) 



dx° A dx° B d 2 k a± d 2 k bl _ je "^ A ^ tJeiTWr J d 2 ka±d 2 h± 
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The value of the factorization scale /i in the parton densities is arbitrary, and since 
from the collinear-factorization standpoint (cf. the Introduction) eq. (|167|) has the same 
accuracy as a LO calculation, the dependence of the jet-production rate ( |178| ) on /i is 



maximal. In eq.(l78) Mueller and Navelet fix the factorization scale at the transverse 
momentum cutoff, /i = p_i_ m in, in order to take the parton densities out of the integration 
over the jet transverse momenta. Then the jet-production rate at fixed x's may be 
directly related to the total parton cross section ( |172| ), 

da 



feff(x%p 2 ±min ) feff(x° B ,p 2 ±m J a gg . (179) 



The total parton cross section at the Born level, in the large- y limit, is given by eq. (|54|) . 
with t integrated over the transverse-momentum cutoff p 2 j_ min - We use it, with the parton 
momentum fractions fl55[ ), in the factorization formula (|179|) to compute the two-jet 



production rate at fixed x's, at the Born level. Normalizing then the two-jet production 
rate ( |179| ), with a gg given by eq. ( |172|) , to the one at the Born level the parton densities 



cancel out. Thus, the if-factor of two-jet production at fixed x's, i.e. the growth rate 
due to the radiative corrections, coincides with the i^-factor of the total parton cross 
section ( |174| ) |2"I |. Then Fig.|9| implies that, could we perform such a measurement and 



were the BFKL approximation correct, we should see a large enhancement in the data 
with respect to the Born-level estimate of two-jet production at fixed x's. 

4.2 Two-jet production at fixed-energy colliders 



In sect. f4.1| we have discussed how, in a hadron collider at variable center-of-mass energy, 



the growth rate of two-jet production at fixed parton momentum fractions may be related 



62 



to the growth rate of the total parton cross section due to the BFKL pomeron. However, 
this measurement may be difficult to implement experimentally, since we do not dispose 
nowadays of a variable-energy collider where such a ramping run experiment may be 
performed^. 

In this section we consider then two-jet production at fixed \fs as a function of the 
jet rapidities ( |177| ). Since at fixed \fs and rapidities, the x's grow linearly with the jet 



transverse momenta (p8|) , the integration over p± in the jet production rate will entail 
a varying contribution from the parton densities, which we cannot avoid. Our goal, 
however, is to examine the parton dynamics and not the parton densities, so we may 



at least fix, or integrate out, the rapidity boost y = (y a + Vb)/^ (cf. sect. [2.3|) , since the 



parton dynamics fl64|) and flBq) does not depend on it. Then we rewrite the factorization 
formula ( |177| ) as 



where the parton cross section is, from eq. ( |167| ) and ( |168| 



rlrr Qn/ 2 00 r°° ( h 2 



^ poo I L™-^ \ 

dk 2 al _dkl ± d(j> ±7rk 3 a± k s b± n f^ Jo \ k\J 

First, we consider the evaluation of the K-factor at fixed yfs, as a function of y. In 
sect.^]5| (cf. Fig.||) the two-jet production cross section da/dydy as a function of y, 
at y = 0, has been computed at the Born level at the energies of the Tevatron and 
LHC colliders. In ref.|M| (cf. also ref. |33| and ||35|| ) the resummation of the radiative 



corrections, using eq.( |180| ) and ( |181| ), has been considered. At Tevatron energies the 

■'A ramping run experiment could be performed in DIS at the HERA electron-proton collider ]5q ], 
by varying the electron-energy loss, but it is not clear whether the rapidity span in jet production in 
DIS at HERA is large enough to show the enhancement effect in the if-factor. 
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corrections turn out not to be appreciably different from the Born calculation, while at 
LHC energies they are substantially larger than the Born estimate, but anyway much 
smaller than what expected from the calculation of the i^-factor at fixed x's. This is due 



to the parton-density contribution to eq. ( |180|) . As we know from the analysis of sect. 2.5 
and Fig.^|, the parton densities fall monotonically as y increases. However, the decrease 
in phase space is more pronounced for the resummed two-jet cross section than for the 
one at the Born level. In order to see it, let us consider the behavior of the factorization 
formulae ( |180| ) and (|60|) as x — > 1. The parton densities have a (1 — x) 13 behavior, with 
/3 positive, as x — > 1. The x's grow linearly with the jet transverse momenta, according 
to eq.fl5"5|) and (^). Trasforming then the integration over the transverse momenta into 
an integration over the x's, we see that the resummed two-jet production ( |18(J| ) goes 
like (1 — x) 2 ^ +1 \ as x — > 1, while the one at the Born level (|60| ) goes like (1 — x) 2/3+1 . 
Accordingly, the corresponding i^-factor falls like (1 — x) as y increases |3 



Another feature of the parton cross section in the BFKL formalism (|167|) is the 
jet-jet decorrelation in transverse momentum p± and azimuthal angle (p. Eq.( |167| ) has 
embodied the correct transverse-momentum conservation, and at the Born level ( |169| ) 
only two partons, balanced in p± and back-to-back in <ft, are produced. However, in the 
large-?/ limit the leading contribution to the eigenvalue ( |156| ) in eq. ( |167| ) comes from the 
n = term (|159|) , for which the partons are uncorrelated in 0; and also the correlation 
in p± fades away as y increases ( |171| ). So eq. (|181|) should allow us to go from one end 



(total correlation at the Born level) to the other (total decorrelation at asymptotically 
large y) and examine how the decorrelation in p± and increases with y. This analysis 



has been performed for the p^-decorrelation in ref. [Q, ||34|| , and for the 0-decorrelation 
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Figure 10: <fi distribution at y/s = 1.8 TeV, normalized to the uncorrelated cross section 
da/dydy. From top to bottom, relative to the peak, the solid lines are the distributions 
at y— 5, 6 and 7. 



in ref . ||34j| , ||35|| , |p7| . From ref.|3J] we report in Fig.|10| the distribution at the Tevatron 
energy as a function of y, at y — 0. The jet transverse momenta are integrated above 
P_i_min = 20 GeV, and the distribution N(<p) is obtained from the two-jet production rate 
da/dydydcf) (|180|) , normalized to the uncorrelated one da/dydy. The renormalization 
and factorization scales are set to fj 2 = k a ±kb±- 



The plot of FigjK] may be understood in terms of the BFKL-ladder picture. As 
the rapidity interval y between the tagging jets is increased, along the gluon ladder 
more and more partons are produced which decorrelate the tagging jets in p± and 0. 
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The decorrelation in p± and as ?/ is increased has been indeed observed in two-jet 
production Hl8| , f58|. However, it is still to be established whether the interpretation of 



this phenomenon in terms of the BFKL ladder is correct, since the decorrelation may be 
also produced via the collinear emission in the parton-density evolution |j8 

4.3 The higher-order corrections 



The jet-jet decorrelation examined in sect.fO raises a few questions. How important 



are in the jet-jet decorrelation configurations where three or more jets are produced? 
On kinematic grounds, we would expect these configurations to be relevant when the 
tagging jets are not correlated. How well are reproduced within the BFKL formalism 
configurations with three or more jets? To answer these questions we must examine the 
kinematics in greater detail. From eq. ( |128| ) we know that in the multiregge kinematics 
transverse- momentum conservation is correctly taken into account, however in the light- 
cone momentum conservation (^) only the leading term is kept. At the Born level, where 
both the produced partons are tagged as jets, we know that this is a good approximation 
for y ^ 3 (cf. sect. 2.5). In this section we examine how good the approximation is when 
we consider higher-order corrections to the Born- level configuration. In order to do that, 
since we do not dispose yet of complete next-to-leading logarithmic corrections to the 



BFKL formalism [59], we use the knowledge of the exact three-parton kinematics and 



matrix elements \TU\. 



In sect. 3.2 we have noted that the exact matrix elements for three-gluon production 



in gluon-gluon scattering ([84] ) reduce in the limit of strong-rapidity ordering (|6~5"D and 
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Figure 11: p± distribution of the forward jet in rapidity, with the transverse momentum 
of the trailing jet fixed at 50 GeV, at y — and at y = 2 and 6. The jet-cone size R cut 
is fixed at 0.7. 



(|66|) to the ones computed through the BFKL formalism (|83|). To examine the accuracy 
of this approximation as well as of the one on the light-cone momentum conservation 
(|68|), we report from ref. |43j the contribution of the three-parton configurations to the 



p± decorrelation. In Fig.|ll| the p± decorrelation is plotted at ^fs = 1.8 TeV as a function 
of the transverse momentum pu_ = k a j_ of the forward jet in rapidity, at a fixed value 
of the transverse momentum p2± = = 50 GeV of the trailing jet in rapidity. The 
solid curves are computed through the large-?/ parton cross section and kinematics 
(|68|), using the factorization formula ( |180|) ; the dotted curves are computed like the solid 
ones, but using instead the exact kinematics ( |63|) for three-parton production; the dashed 



67 



curves are computed through the exact matrix elements |TI| and kinematics (j63|) , using 
the three-parton contribution to the exact factorization formula ( |1 76|) . All the curves 
show an infrared divergence at k a j_ ~ ki,±, where the third parton becomes soft. In 
the exact calculation also collinear divergences are present, when two of the final-state 
partons become collinear. These are disposed of by discarding configurations^] where the 

1 /2 

distance R = [{yi — yj) 2 + (4>i — 4>j) 2 } between two partons on the lego plot is smaller 
than the jet- cone size R cu t- 



From Fig.[Tl] we see that the error in using the large-?/ approximation grows with the 
imbalance in transverse momentum of the tagging jets. The dotted curves, which are 
theoretically inconsistent, are plotted just to show that while at small y's the error is 
distributed between the approximation on the matrix elements and the one on the parton 
densities, at large y's most of the error comes in using the large-;?/ approximation ( p8j ) 
in the parton densities. To understand it, we recall that in the three-gluon cross section 



(pq) the intermediate gluon ki is produced with equal probability over the rapidity range 
determined by the tagging jets. This is also true at the hadron level if we neglect the 
contribution of the intermediate gluon to the light-cone momentum conservation (|68|) . 
However, we see from the exact kinematics fl63|) that this can be a bad approximation 
if ki± is not small, particularly when the intermediate gluon is close in rapidity to the 
tagging jets. In this case using the exact kinematics in the parton densities produces a 
large suppression, so the rapidity range that the intermediate gluon may effectively span 

is reduced substantially. The BFKL (solid) curve neglects this effect, and so it grossly 

fc In a complete 0(a\) calculation these configurations would be counted as two-jet events with k a ± = 
kb±, however they may be neglected in our analysis, since we are only interested in the modifications 
induced by the third parton to the large- y kinematics when k a j_ =/= kb±. 
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overestimates the cross section^. However, at k a ± ~ kb± the transverse momentum of the 
intermediate parton is small, so its contribution to the x's in (|63|) can be safely neglected. 



This entails that the BFKL approximation works at its best when the tagging jets are 
balanced in p±_ and back-to-back in <ft, and may have large errors, even at large y's, when 
the jets are not correlated. This picture suggests that the tails in the distributions 
of Fig. [Ty may be overestimated, in other words that the decorrelation given by the 
BFKL approximation is too strong. This seems to be confirmed by the preliminary data 
of the DO Collaboration 



A remedy to this problem is to use the knowledge of the exact three-parton config- 
urations and the ambiguity in the definition of rapidity in the BFKL approximation. 
We identify the rapidity Ubfkl = ln(s/£r[) in the BFKL ladder, with k\ a typical 
transverse-momentum scale, with the kinematic rapidity y = y a — y^. However, the 
BFKL formalism being a leading logarithmic approximation, ysFKL is defined only up 
to a factor which is subleading at large rapidities. In ref. ^| this ambiguity has been 
exploited in order to define an effective rapidity y which differs from the kinematic one y 
by subleading terms, and which keeps into account the exact three-parton configurations 
described above. Using y, instead of y, in the BFKL resummation for the <fi decorrelation 



has yielded a much better agreement with the data |p8 |. However, this is only a 
phenomenological improvement, since from the theoretical point of view to use y or y in 
the BFKL resummation is equally valid. 



'We have found that even at LHC energies this discrepancy is not negligible. For example, for the 
configuration of Fig.|ll| with y = 10, k a± = 20 GeV, k b± = 50 GeV, R cut = 0.7, the BFKL curve 
overestimates the exact calculation by almost a factor 2. 
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5 Conclusions 



In these lectures we have analysed the BFKL formalism, which resums the leading log- 
arithmic contributions, in m(s/i), to the radiative corrections to two-parton processes, 
and we have applied it to inclusive two-jet production at large rapidity intervals. 

After sketching the parton-model and factorization ideas in the Introduction, all the 
ingredients necessary to cook up two-jet production at leading order in a s are listed in 
sect. 2; then the two-parton dynamics is analysed in detail in the large-y limit (sect. 2.4) 
and it is shown that the large-?/ approximation to two-jet production agrees well with 
the exact calculation for y ^ 3. 

In sect. 3 the tree- level three-gluon amplitudes ( ff9|) and the one-loop corrections to 
the two-gluon amplitudes ( |102j ) are computed in the large-?/ limit, and an ansatz is made 
on the elastic ( |106[ ) and the multigluon ( |108| ) amplitudes in the leading logarithmic ap- 



proximation; its self-consistency is shown by using the multigluon amplitudes to evaluate 
the elastic amplitude for gluon-gluon scattering ( p.06| ). In evaluating the elastic ampli- 
tude, the BFKL equation ( |136| ) which describes the evolution in transverse momentum 
of the gluon ladder exchanged in the t channel is introduced. The solution of the equa- 
tion for color-singlet exchange at t — ( |152| ), corresponding to the exchange of a (cut) 
BFKL pomeron, is then related to the total cross section for gluon-gluon scattering in 



the large- y limit (|166|) . The analytic structure of the total cross section is examined in 



sect. 3.8, and its growth rate with s, due to the BFKL pomeron, is stressed in Fig. 9. 
In sect. 4 the BFKL formalism is applied to estimate two-jet inclusive production at 
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large y's, starting with Mueller-Navelet proposal which relates the growth rate of two-jet 
production at fixed x's ( |1 79| ) to the one of the total parton cross section; in sect. 4.2 
then two-jet production at fixed-energy colliders at large y's is considered, and the jet- 
jet decorrelation in transverse momentum as a function of y is proposed as a signature 
of BFKL dynamics. Finally, in sect. 4.3 knowledge of the exact next-to-leading order 
kinematics is used to improve the phenomenological analysis of the jet-jet decorrelation. 

Even though the parton dynamics described by the BFKL formalism is rather so- 
phisticated, a few caveats are in order in dealing with its applications to jet production: 
the formalism has, from the collinear-factorization standpoint (cf. sect. 1), the same 
valence as a leading-order calculation, accordingly it has the maximal dependence on 
the renormalization and factorization scales; the rapidity interval to use in the BFKL 
resummation is defined only up to subleading terms, which may be quite important at 
the energies of the Tevatron and LHC colliders (cf. sect J4 .31); the jet structure, being a 
subleading feature of jet production [)i"3|| , cannot be resolved by the BFKL formalism. 
These considerations and the results of sect. 4.3 seem to suggest that a complete calcu- 
lation of next-to-leading logarithmic corrections to the BFKL formalism should improve 
substantially the analysis of jet physics at large rapidities. 

Acknowledgements 

I am very grateful to James Bjorken for triggering and nurturing my interest in this 
topics, to Michael Peskin with whom much of sect. 3 was worked out, to Carl Schmidt 
for a fruitful collaboration. I also wish to thank Joachim Bartels, Lev Lipatov, Alfred 



71 



Mueller, Mark Wiisthoff, Peter Zerwas for many useful discussions, and Enrico Predazzi 
for encouraging me to write this whole thing up. 



A More on the BFKL equation at t = 



Here we consider the BFKL equation at t — ( |146| ), without making the replacement 
( |1 13| ) in the virtual radiative-corrections term. Substituting the solution ( |152| ) into the 
homogeneous equation associated to eq.( |146| ), we obtain the spectrum of the eigenvalues 



co>(V, n) 



a«N r . 
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2 Re dx 

o 1 



n|-l 
X 2 



X 



dx 

o 1 — x 



-I dx — 

2 Jo x 



with x defined as in eq.( |154j ). The eigenvalue has divergences at x = 
x = 1 <-> (q 2 — qi)\ = 0. In order to regulate them, we require that 



182) 



12± = 0) an( l 



Q2± > ^ 



2 I 2 

x > n /q 



i± j 



183) 



(l2- qi)±> fJ 2 x<l-/j,/q 



i± 



with yU 2 an infrared regulator. Then in the last term of eq. ( |182|) we move the singularity 
from ql L = to (q 2 - qi)± = 0, 



"1-^/92 1 

, dx — 

2 J^/ql X 



1-^/92 



dx ■ 



1 — x 



V /? 2 

and we obtain a spectrum which is well behaving in the integration limits 

r 1 , x^ +lv - 1 

/ dx 

Jo 



uj(v, n) 



2 Re / dx- 



184) 



185) 



7T JO 1 — X 

and agrees with eq.( |156| ). The calculation we have outlined here is equivalent to the one 
done in sect. |3T7| indeed the regulation of the infrared divergences that we do in eq.( |183| ) 
and ( |184| ) is done in sect. [37^ by eq. ( |143j ). 
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B The eigenvalue uj(y,n) 



From the properties of the ip function the largest values of the eigenvalue in eq.( |156| ) 
are at \n\ — 0. In particular for v = we have 



u)(v = 0,n) = -4- 



, a s N r , 
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m j 

— In 2 + for Inl = 2m . 

1 m 1 



186) 



2 fc=i fc 



for |n| = 2m + 1 , 

So we fix |n| = and look at the dependence of uj(u,n) on v. Using the expansion []52 



V»(l + y)-^(l) = J2(-l) k ak)y k 

k=2 



fc-1 



187) 



valid for small y, with C{k) the Riemann ^-function, and the doubling formula, 

2^{2y) = 2m2 + ^(2/)+v(2/+^) , (188) 



we obtain 



u>{v, n = 0) = 2- 



7T 



2 In 2 + ^(-l) fc (2 2k+l - 1) ((2k + 1) z/ 2fc 



k=l 



189) 



C Inserting pomeron amplitudes on Born diagrams 



We show here a simple way of inserting one or more BFKL ladders on cut diagrams at 
the Born level. Let us consider the Born cross section for the production of gluons k a 
and k b , in the large-?/ limit (|54D, rewritten as 



' 99 



L/ A a s r2/ 



d?k a i_d?k]ji_ 2kaj_ k 2 j_ 
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(190) 



Then we insert the (cut) pomeron off-shell amplitude f(k a ,kb,y) ( |168| ) on the Born 
diagram, and assume that the cross section for the production of two gluons, resummed 
to all orders of a s y has the form, 

The coefficient C here is unknown, and is determined by requiring that at the Born level 
eq. (|191|) reproduces eq. ( |190|) . At the Born level there are no radiative corrections, so the 



singlet solution ( |152| ) is given by the inhomogeneous term of eq.( [f47| ), 

fi(k a , h) = $ 2 (k a ± + h±) , (192) 



where as usual k a ± = q al _ and k b ± = —qb±. Substituting eq.( |168| ) and (|192| ) into eq. ( |191|) 



and comparing with the Born cross section ( |190| ), we obtain C = C\, and eq.( |191| ) agrees 
with the resummed parton cross section (|167|) . 

In the same way we may derive the BFKL corrections to the cross sections for pro- 
duction of two heavy-quark pairs in photon-photon scattering fl2"7| , for forward-jet pro- 
duction in DIS |p6[ , for (Higgs-boson + one-jet) production in gluon-gluon scattering 
provided the respective Born cross sections are known in the large-?/ limit. 



Now suppose that we want to derive the cross section for the production of three 
gluons k a , k b and k c , with rapidity ordering y a ^> y c ^> y&, assuming that between each 



pair of gluons (k a , k c ) and (fc c , kb) a BFKL ladder is exchanged |p0[ . The corresponding 
Born cross section may be derived from eq.(|88|), or from eq. (|191| ) with the singlet solution 
( |152|) corresponding to the emission of a real gluon in the rapidity interval y a — yb (cf. 
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Figure 12: Three-gluon production with exchange of two BFKL ladders, represented 
schematically by the double-lined squares. 



the kernel of eq. (|147|) ), 



fi{k a ,h) 



1 



1 



a s C A 



tt 2 (k a -k b y 2(i-iy 



and we obtain 



da, 



99^999 



cWs 



■5 2 {k a ± + h± + k 



c± 



(193) 



(194) 



d 2 k a± d 2 k b± d 2 k c± dy c ^ 2 k 2 a± kl ± k 2 c ^ 
where the rapidities of gluons k a and k b are as usual fixed by the parton momentum frac- 



tions. Guided by eq.([194j), we assume that the cross section for three-gluon production 
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with the insertion of two BFKL ladders has the form (Fig.[T2|), 

dfigg->ggg (lg5) 



d 2 k a± d 2 k b± d 2 k c± dy c 

>,3 



ex f 

D 1,2 A t.2 / d 2 Qcd 2 q d 5 2 (q cl _ - q dl _ - k c± ) f (q a , q c , y ac ) f (q d , q b , y cb ) , 

K a± K b± K c± J 

with y ac = y a — y c and y cb = y c — yb- In order to determine the coefficient D we take 
the singlet solution at the Born level ( |192| ), and insert it into the pomeron amplitudes 

ai qa Vac 

) and f (q d ,q b ,y c b)- By comparing then the result with the corresponding 
Born cross section ( |194j) , we obtain D = 2C\/tt 2 . Thus, eq.( |195| ) agrees with the three- 
gluon production cross section with exchange of two BFKL ladders computed in ref. 
||60|| . This procedure is straightforwardly generalizable to the production cross section 
for n rapidity-ordered gluons, with exchange of a BFKL ladder between each pair of 
nearest-neighbor gluons. 

D The diffusion equation 

In order to understand the asymptotic behavior ( |171| ) of the parton cross section, which 
follows from taking the leading singularity ( [L60| ) of the partial-wave solution in the 
complex plane of /, we must consider the recursive relation (|135|) at a very high value 
of the order parameter n in the iteration procedure. Accordingly, we take the BFKL 
equation ( |147|) at a high value of n, for which the inhomogeneous term is not important. 



Then from eq.(148) we obtain 



(l-l)t\ qi ,k) = (196) 

a„N r . 



7T 2 



/ d 2 q 2± 



ft 1] (Q2, k) - (q 2 Jq 2 ± ) ft l \^ k) + QUt^iQi, k) 



(?i - 12)1 ql± [ql±_ + (qi ~ ^2) 
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Next, we take the solution ( |152|) , averaged over the azimuthal angle, and consider only 
its dependence on qi±, 



fi n \qi) ~ (ql±r 1/2 M^i) , with Ai = Hqfjm 2 ) 



(197) 



and we replace it in eq. (|196|) , in the limit (n| ln(g 2 _|_/g2_i_)|) «C 1, for which the produced 
partons are approximately balanced in q±. Then we expand the solution ( |197|) taking 
the configuration where the partons are balanced in q± as a stable point, 

(A 2 -A 1 ) 2 9 2 ^_i(A 1 ) 



tpn-l(h) = ^n-l(Al) + 

and for the recursive relation ( |196| ) we obtain 

' e (A a -A 1 )/2 _ i 



(198) 



(Z-l)^(Ai) 



a s N c rx 



7T JO 



II - e A2 ^ Al 



+ 



1 



Vl + 4e 2 ( A2 - Al ) 



^n-l(Al) 



+ 



(A 2 -A x ) 2 e ( A2 ~ Al )/ 2 #Vn-l(Al) 



(199) 



2 |l-e A 2- A i| 8\\ 
The integration over the first two terms on the right-hand side yields the eigenvalue ( p.56|) , 
of which we take the asymptotic solution ( |159| ). Neglecting then subleading terms, the 
recursion relation becomes 



(I - l)^(Ai) = A^n-i(Ai) H / dz- 



)V> (200) 
7r Jo e* — 1 aAf 

with A given in eq. (|161|) . Using the integral representation of the Riemannn ^-function 
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C(*) 



(2* - i)r(x) Jo 



e 2 *-l 



(201) 



and B as given in eq.( 161 ), and going to continuous values for n, we obtain a diffusion 
equation, 

(1-1) ^ = \A-{1-1)] iin, A) + B ^> . (202) 



dn 
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This is a Schrodinger-type equation, with n playing the role of time, and with a constant 
potential, A — (I — 1). Therefore the solution is separable U6lj , ip(n, A) = f(n)g(X). At 
the singularity I — 1 + A the diffusion equation is, 



gV<n, A) = <9 2 ?/;(n,A) 



(203) 



0n OX 2 ' 

with diffusion rate A ~ n 1 ^ 2 . Supposing that at the time n = the wave function has a 
Gaussian shape, 



ip(n = 0, A) 



1 



exp 



A 2 



2a 2 



(204) 



we obtain the solution as a wave packet spreading in time |4"6 |, in qualitative agreement 
with eq.( |171| ) if we make the correspondence Ay ~ n, 



(205) 



where we have neglected the initial width, Ao 2 <ti 2Bn. 
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